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^ SUMMARY 

A detailed  formulation  of  the  equations  of  motion  of  a deformable  aircraft 
is  given.  The  development  is  from  Lagrange’s  equations  for  an  inertial  frame, 
and  is  made  in  terms  of  the  position,  orientation,  force  and  inertia  properties 
of  narrow  strips  of  the  aircraft  which  lie  fore  and  aft  in  the  unperturbed  state. 
The  latter  is  one  of  constant  linear  velocity  and  zero  angular  velocity. 

Particular  account  is  taken  of  the  deformation  and  loading  in  the  unperturbed 
state. 
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1 INTRODUCTION 

There  are  two  ways  of  trying  to  understand  a complicated  problem.  One  is, 
as  it  were,  to  stand  back  and  see  the  thing  as  a whole  without  being  confused  by 
all  the  detail  - to  look  at  the  wood  and  not  the  trees.  The  other  is  to  pick 
on  a basic  unit  - the  tree  of  ^he  phrase-  and  look  at  that  first.  In  each  case 
one  has  ultimately  of  course  to  bring  in  the  details  - to  look  closer  at  the  wood 
or  to  build  up  the  wood  from  the  trees.  In  Ref  1 the  first  approach  was  adopted; 
but  a recent  note  by  Baldock^  showed  that  there  are  some,  and  they  belong  to  the 
category  of  those  who  do  it  rather  than  those  who  show  you  how  to  do  it,  who  want 
to  use  the  other  route.  The  consequence  in  this  paper,  which  is  a considerable 
generalisation  of  what  Baldock  did,  written  in  a notation  which  harmonises  with 
that  of  Ref  1 (and  indeed  of  Refs  3 and  4) . 

2 THE  DATUM  MOTION 

We  wish  to  consider  perturbations  from  a datum  motion.  What  then  should  be 
the  unperturbed  state?  How  general  a motion  should  it  be?  What  restrictions 
should  be  placed  upon  it?  One  view  would  be  to  select  it  to  suit  flight  conditions 
encountered  in  practice;  but  alternatively  one  could  choose  it  to  ensure  the 
minimum  of  complication  in  the  analysis,  and  then,  if  necessary,  generalise  when 
required.  We  have  taken  the  latter  approach  and  in  particular  have  said  that  the 
datum  motion  shall  be  such  that  the  shape  of  the  aircraft  remains  constant 
throughout.  In  particular,  this  ensures  considerable  simplification  of  the  aero- 
dynamics and  kinematics  of  the  system.  To  ensure,  in  general,  that  there  is  no 
change  of  shape  during  the  datum  motion  we  therefore  specify  that: 

(i)  the  aircraft’s  mass  and  mass  distribution  are  constant, 

(ii)  the  atmosphere  and  the  earth's  gravitational  field  are  uniform,  and 

(iii)  the  motion  is  one  of  constant  linear  velocity  and  zero  angular 
velocity  relative  to  earth  fixed  axes. 

We  also  assume  that  the  aircraft  is  symmetric  in  the  unperturbed  state. 

Consequently,  we  take  a set  of  constant-velocity  axes  °cxcyczc  which  has 

zero  angular  velocity  and  constant  linear  velocity  in  a direction  in  the 

plane  0 x z , and  is  such  that  it  is  fixed  in  the  aircraft,  0 x being  fore 
c c c c c 

and  aft,  throughout  the  datum  motion.  Thus  the  datum  motion  is  defined  by: 


where  0^  is  the  angle  of  inclination,  y ^ is  the  angle  of  climb,  the  angle  of 
bank  is  zero;  and,  with  no  loss  of  generality  for  our  purpose,  we  have  taken  the 
nose-azimuth  angle  and  the  angle  of  track  both  also  to  be  zero*.  Thus  the 
constant-velocity  axes  coincide,  during  the  datum  motion,  with  the  body-fixed 
axes  used  to  define  the  attitude  of  the  aircraft.  One  could  choose  the  latter 
axes  to  be  the  principal  axes  of  inertia  of  the  aircraft  and  consequently  obtain 
some  simplification  in  the  expressions  for  the  reversed  effective  and  generalised 
gravitational  forces  at  the  expense  of  having  to  find  out  where  the  principal 
axes  of  inertia  are  (of  Appendix  B) . Another  choice,  which  may  appeal  to  some, 
would  be  to  take  body-fixed  axes  which,  for  the  particular  datum  motion  being 
considered,  coincided  with  the  body-path  axes  having  the  same  origin.  In  this 
case  one  would  have  the  angles  of  climb  and  inclination  equal  (y^  ■ 0^) . 

3 THE  SEMI-RIGID  MODEL 

We  assume  the  aircraft  consists  of  a number  of  strips  each  of  which  is 
based  on  a mean  line  which,  in  the  datum  motion,  is  normal  to  °cyc  ♦ passes 

through  the  strip  centre  of  gravity,  and  is  more  or  less  fore  and  aft.  Each 

strip  is  made  up  of  two  rigid  portions**  joined  by  a hinge  at  a point  on  the  mean 

line.  The  plane  separating  a strip  from  an  adjacent  strip  is  not  necessarily 

normal  to  0cyc  , but  instead  is  chosen  to  be  roughly  normal  to  the  external 
surface  of  the  aircraft  at  their  line  of  intersection.  We  also  arrange  that  the 


* The  three  attitude  angles  (angle  of  bank,  angle  of  inclination,  nose-azimuth 
angle)  define  the  attitude  of  the  aircraft,  ie  of  some  axes  fixed  in  the 
aircraft,  relative  to  normal  earth-fixed  axes.  The  two  flight-path  angles 
(angle  of  climb,  angle  of  track)  define  the  direction  of  flight  relative  to 
the  same  axes  (of  Ref  3,  sections  5 and  6). 

**This  image  includes  the  case  of  a single  rigid  portion  as  the  case  where  one  of 
the  two  portions  is  of  zero  length. 
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mean  line  of  a strip  does  not  extrude  from  the  strip  through  either  of  its 
separating  planes.  Such  a strip  models  a wing  (or  tailplane  or  fin)  chordwise 
strip  with  a control  surface;  or,  with  rather  more  engineer's  licence,  even  an 
engine  pod  or  a store  or  a fuselage. 

In  Ref  1 an  intermediate  frame  of  reference,  defined  by  the  no-deformation- 
body-fixed  axes  ®nxnynzn  » was  u8e<^  in  the  specification  of  the  perturbed 
position  and  shape  of  the  aircraft.  In  the  present  context  it  is  convenient  also 
to  use  sets  of  strip-fixed  axes,  one  set  being  associated  with  each  strip.  Let 
the  strip-fixed  axes  for  a typical  strip  be  0 .x  .y  .z  . where  0 .x  . is  the 
mean  line  mentioned  above,  the  hinge  is  on  0 .x  . , ie  at,  say, 

S X S X 


(3-1) 


and  the  axes  are  fixed  in  one  of  the  rigid  portions  of  the  strip  such  that  the 
hinge  axis  of  rotation  is  parallel  to  the  direction  » an<*  ®sixsiZsi 

a mean  plane  between  the  planes  dividing  the  strip  from  the  adjacent  strips. 

Thus  for  a wing  section  with  a trailing  edge  control  one  would  take  0g^  to  be 
some  reference  point  such  as  the  quarter  chord  point,  °s£xs£  to  be  forward 
along  the  chord  line,  and  x^g^  to  be  the  (negative)  value  of  xg^  at  the 
control  surface  hinge*.  We  will  call  the  portion  of  a strip,  in  which  the  strip- 
fixed  axes  are  fixed,  the  main  part;  and  the  other  portion,  when  it  exists,  the 
flap  ( cf  Appendix  B) . 


As  in  Ref  1 we  represent  the  perturbation  of  the  aircraft,  from  the 
unperturbed  position  and  shape  it  would  have  had  at  the  same  instant  dur’ *r  the 
datum  motion,  as  being  made  up  of  rigid  body  translations  and  rotations  which 
move  a set  of  body-fixed  axes  from  coincidence  with  the  constant-velocity  axes  to 
coincidence  with  the  no-deformation-body-fixed  axes,  followed  by  some  further 
perturbations  which  we  call  deformations.  The  first  part  of  these  perturbations 
is  therefore  (of  Ref  1)  made  up  of  the  two  successive  steps: 

Co)  (c)  ” (c)  • • • 

(i)  translations,  as  a rigid  body,  xj  , y|  , z|  in  the  directions 
of  the  respective  constant-velocity  axes; 

(ii)  successive  rotations  p,  6,  $ about  the  carried  axes  Oz,  Oy,  Oz 
where  Oxyz  are  the  above-mentioned  body-fixed  axes. 


* The  control  surface  hinge  and  the  hinge  of  the  strip  model  do  not  necessarily 
coincide.  If  the  control  surface  has  a swept  hinge  and  the  strip  is  fore  and 
aft  they  will  not. 
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The  second  part  of  the  perturbations  - the  deformations  - we  write  in  a somewhat 
different  form  from  that  of  Ref  1. 

The  transformation  from  the  no-deformation-body-fixed  axes  to  the  strip- 
fixed  axes  for  the  ith  strip  is  achieved  by: 

(a)  Translations  xfn^ , yfn\  zfn^  in  the  directions  of  the  respective 
no-deformation-body-fixed  axes  where  these  translations  are  compounded  of  the 
datum  state  separation  between  the  strip  reference  point  0 . and  the  point  0 
and  a deformational  contribution,  viz: 


f (nJl 

x; 

l 

Xif 

+ K 

ql 

(n) 

yi 

yif 

• 

e 

Lzi  J 

-if 

_qn_ 

K is  a modal  matrix  whose  elements  are  functions  of  the  strip  being  considered, 
and  qj  -+■  qn  are  the  generalised  coordinates  for  the  deformational  freedoms. 

For  example  if,  for  a wing,  the  strip  reference  points  are  at  the  quarter  chord 
point,  then  K describes  the  deformational  shapes  of  the  quarter  chord  line  in 
the  various  modes*. 


(b)  Successive  rotations  i|^,  6^,  <{k,  about  the  carried  axes,  which  are 
written  in  a form  similar  to  (3-2) : 


h 

- 

"♦if 

+ F 

qf 

0. 

l 

eif 

* 

*i 

l—  -J 

♦if 

L. 

_qn_ 

(3-3) 


Thus  F is  a modal  matrix.  It  can  be  thought  of  as  a 'torsional'  modal  matrix 
along  with  K as  a 'flexural'  modal  matrix.  The  condition  that  a strip  mean  line 
°sixsi  norma^  to  °cyc  during  t*'e  datum  motion  is  satisfied  by  making 
= 0 , and  this  does  not  impose  any  other  restriction,  as  we  will  show 
following  equation  (3-21).  For,  say,  a wing  whose  reference  axis  (line  of  strip 
reference  points)  is  parallel  to  the  plane  Onxnyn  (ie  Z£f  • constant),  whose 

* All  the  values  of  K can  be  thought  of  as  a three-dimensional  array  of  numbers. 
It  is  this  array  which  describes,  in  this  case,  the  deformational  shapes  of  the 
quarter  chord  line. 
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only  flexibility  is  in  torsion  about  the  reference  axis  ( \ie  the  matrix  K is 

zero  for  all  wing  strips),  and  whose  strips  are  fore  and  aft  < [ie  0 .x  . parallel 

si  si 

to  °nxnzn)  during  the  datum  motion,  and  will  be  zero  and  6^f  will 

be  the  sum  of  the  strip  jig  'incidence'  and  the  strip  'incidence'  due  to  twist  in 
the  datum  motion*. 

In  addition  to  the  displacements  and  deformations  represented  by  these  two 
axes  transformations  we  have  deformation  in  the  strip  itself.  The  main  part  of 
the  strip,  by  definition,  cannot  be  displaced  relative  to  the  strip-fixed  axes, 
but  the  flap  part  can.  We  therefore  have  on  the  flap  part  of  the  strip** 


si 

^si 

T 

I + (6.  - 6.,) 
l if 

"0 

0 

f 

"x  . - x,  r 

sie  hsi 

y(s) 

J SI 

0 

< 

0 

0 

0 

► 

y . 

•'sie 

Z(S) 

_ Sl_ 

_ 0 _ 

• 

0 

0_ 

4 

z • 

sie  -1 

(3-4) 


where  6^  is  the  angle  of  flap  rotation,  6.^  is  the  flap  rotation  in  the  datum 

motion,  x . , etc  are  the  values  of  x . , etc,  when  there  is  no  flap  rotation, 
S 16  S X 

and  is  t^>e  axes  transformation  matrix  (attitude  deviation  matrix)  for  the 

single  rotation  6^  about  0s.y  . which  is  given  byt 


cos  6.,  0 

if 

0 1 

sin  6.,  0 

if 


- sin  6., 
if 


cos  6.. 
lf 


(3-5) 


* Strictly  speaking  these  angles  should  be  called  angles  of  inclination  (in  the 
no-deformation-body-fixed  axes  reference  frame)  rather  than  angles  of 
incidence  ( of  Ref  3,  sections  5 and  6). 

**  We  have  used  the  general  symbol  6 , specified  by  Hopkin  (Ref  3,  section  8) 
for  a motivator  rotation  rather  than  decide  the  primary  purpose  of  our  flap 
(say  to  produce  a rolling,  pitching  or  yawing  moment)  and  take  the  appropriate 
symbol  (£,  n or  C) . 


* Thus  d6.f  (P<sif) 

if 


0 0 
0 0 
1 0 


j I 
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(For  a general  rotation  the  axes  transformation  matrix  S is  of  the  form 
S = R.PqY,  - of  Ref  4,  Appendix  A.)  The  above  expression  (3-4),  assumes  that 

<f>  0 \jJ 

6^  - 6^  is  small.  Putting  <5^  = <5^  , as  in  the  datum  motion,  we  therefore 
find  that*,  on  the  flap. 


"(us)" 

Xsif 

s 

^Sisi 

T 

+ P* 

6if 

Xsie  xhsi 

(us) 

ysif 

0 

ysie 

(us) 
x sif  _ 

0 

z . 

_ sie 

(3-6) 


The  perturbation  in  the  flap  rotation  is  related  to  the  generalised  deformational 
coordinates  by  a flap  modal  vector  f : 


6.  - 

l 


(3-7) 


and  so,  on  the  flap 


r (s  5i 

X . 

SI 

f^s 

” (usf 
xsif 

T 

0 

0 

1 

X . - X,  . 

sie  nsi 

fT 

qi 

(s) 

y • 

■'si 

(us) 

ysif 

0 

0 

0 

y . 

J sie 

• 

• 

Z(s> 

- S1_ 

(us) 
_ sif  _ 

-1 

0 

0 

z . 
sie 

_qn 

Combining  the  successive  deformations  (equations  (3-2),  (3-3)  and  (3-8))  we 
find  that,  for  a point  on  the  flap,  its  position  relative  to  the  origin  of  the 
no-deformation-body-fixed  axes  is  given  by 


* The  superscript  (us)  is  here  introduced  to  show  that  these  are  the  coordinates 
referred  to  the  unperturbed  (datum  motion)  orientation  of  the  3trip-fixed  axes. 


1 
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(n) 
x . 
m 
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(n) 

# 

y . 
7ni 

yif 

• 

ni_ 

_Zif_ 

_qn 

__  __ 

— 

— 

— — 

x0^ 

S3  f 

T 

+ K 

6if 

0 

0 

1 
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sie 
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(3-9) 


where  . (=  R^.PgjY^)  is  the  axes  transformation  matrix  for  the  change  in 

orientation  from  the  no-deformation-body-fixed  axes  to  the  strip-fixed  axes 
{of  Ref  4,  Appendix  A) 


and 


Writing 


1 0 0 
0 cos  0^  sin  <(>j 

0 - sin  4>£  cos  <f>j 


cos  6 ^ 0 - sin  0^ 

0 1 0 
sin  0^  0 cos  0^ 


cos  sin 

- sin  cos  iJk 

0 0 


0 

0 
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a7 

l 

ei 

"♦i" 
0.  - 

i 

hi 

0if 

= Q*  F 
♦if 

~ql" 

_Yi 

.♦i- 

*i£ 

_qn 

(3-10) 


(3-11) 


(3-12) 
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where 


1 0 

0 cos  <(>.,, 
if 

0 - sin  <)Kf 


sin  0 . , 
if 


(3-14) 


sin  <f . , cos  0 . 
if  if 

cos  A...  cos  0... 
if  if 


it  can  be  shown  that* 


>i  * ('  - *0‘ 


where 

A 

n 

0 

■yi 

h 

• 

1 

Yi 

0 

- a . 

l 

a. 

i 

0 

(3-15) 


(3-16) 


Thus  to  first  order,  making  use  of  the  fact  that,  for  matrices  of  the  latter 
type. 


\rxi = _Axp 


(3-17) 


we  have 

r»i 


x.,  + S' 

if  4 


(us) 

ysif 


+ K - A . \Q.  F + S*  P*  0 0 1 Tx  . - x,  .]  fT  |q; 

♦if  x0^  ♦if  ♦if  5if  Sie  ^S1  1 

. Slf  0 0 0 V • 

ysie  i 

- 1 0 0 z . q 


(3-18) 


* Thus  the  angles  ct£  . . . are  an  approximation  to  the  Euler  angles  which  produce 
the  orientation  transformation  from  the  unperturbed  position  of  the  strip-fixed 
axes  to  their  perturbed  position  when  the  perturbations  if,  0,  if  are  zero 
(of  equation  (4-40)). 
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This  expression  gives  the  position  of  a particle  on  the  flap  part  of  a strip 
relative  to  the  origin  of  the  no-deformation-body-f ixed  axes  and  resolved  along 
th  se  axes.  For  a particle  on  the  main  part  of  a strip  one  merely  takes  the 
modal  vector  f to  be  zero.  This  equation  compares  with  expression 


r*(n)i 

n 

— — 

“ ““ 

= 

xf 

+ R 

ql 

(n) 

y 

'ti 

yf 

• 

(n) 

z 

n 

_zl 

_qn 

used*  in  Ref  1 (equation  (1)).  Thus  for  the  strip  model,  the  modal  matrix  R 
has  the  form: 


R 


T T T 

K - S*  A , . Q,  F + ST  Pt 

<p  . c (us)  <p  . c . , <5 . 

if  x .,  if  if  if 

o 

0 

1  

v — V 

"tie  hsi 

Slf 

0 0 0 

y . 

^ sie 

- 1 0 0 

z . 

- _ 

sie  _ 

(3-20) 


at  the  ith  strip,  and  the  unperturbed  Coordinates  of  a particle,  relative  to 
and  resolved  along  the  body-fixed  axes,  in  the  datum  state,  are: 


xf 

* 

xif 

T 

+ 

♦if 

(us) 

xsif 

yf 

yif 

(us) 

ysif 

_zf_ 

_zif_ 

(us) 

_ sif  _ 

for  a particle  on  the  ith  strip.  If  this  strip  is,  as  required,  in  a plane 

normal  to  0 y then,  for  a point  on  0 .x  . (ie  y .,  = z . = 0) , v,  must 
c-'c  v si  si  if  s if 

be  constant.  Consequently  the  12  term  of  must  be  zero  and  so 

(of  equations  (3-10),  (3-11)  and  (3-12)) 


cos  sin  = 0 


(3-22) 


* There  is  however  a difference  of  some  significance:  (3-19)  is  a precise  state- 
ment of  the  deformation  used  in  Ref  1 while  (3-18),  for  our  present  model,  is 
only  a first  order  approximation  in  the  generalised  coordinates  q. 

(of  Appendix  A) . 


I 
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This  is  satisfied  by  0^  = (n/2)  or  = 0 . However,  any  transformation 

with  0^  = (tt/2)  is  equivalent  to  another  transformation  with  = 0 , and 

so  we  will  always  take 


This  makes 


\ Pfl 

*if  6if 


*if 

= 0 

(3-23) 

cos 

6if 

0 

- sin  6,r  H 
if 

• 

sin  <p . 
if 

sin 

0if 

cos  $if 

sin  p . r cos  6 . .. 
if  if 

cos  <(>if 

sin 

eif 

- sin  6 * r 
if 

cos  cos  0., 

if  if 

(3-24) 


There  is  one  other  expression  for  the  position  of  a particle  on  a typical 
(the  ith)  strip  that  we  may  require,  and  that  is  for  its  position  relative  to 
the  origin  of  the  constant-velocity  axes  and  resolved  along  those  axes.  This  is 
easily  seen  to  be 


|>>1 

c 

xif 

„T 

+ V 

vif 

>*>■ 

Slf 

y(c) 

J c 

yif 

(us) 

ysif 

z(c) 

C 1 

zif 

(us) 

Z • r 
Slf 

+ R 


X * 
— \ 
o 

" (\f  + S*if\(us)Shf) 

p 

(c) 

yl 

0 

Z(C) 

Z1 

P 

(3-25) 


where  R is  given  above  (equation  (3-20)). 

In  this  section,  we  have  introduced  the  generalised  coordinates  q 


/ \ / \ / \ In 

To  complete  the  set  we  take  |xj  y^  zj1  $ 0 i|>}  as  the  other  six  generalised 


coordinates  Iq  . 

1 n+1 

capriciously. 


q g}  and  subsequently  we  will  use  either  notation 


1 


II  Mr  ill' in"  IT-'*- 
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4 THE  AERODYNAMIC  FORCES 

4. 1 The  velocities 

The  main  reason  for  the  adoption  of  the  strip  model  described  in  section  3 
is  of  course  to  facilitate  the  use  of  strip  theory  aerodynamics  so  beloved  by 
many.  The  air  forces  on  a strip  will  depend  on  the  velocity  of  the  air  relative 
to  the  strip.  The  velocity  of  a particle  referred  to  the  strip-fixed  axes  is 


"Ol 

u . 

• (s) 
x . 

+ 

7s?  - ssV  ) 

r<s>i 

X . 

+ 

mi 

SI 

♦ . 

V 4>-  J 

SI 

1 

1 

\ 1 1 / 

(s) 

• (s) 

V 

4 

y(s) 

•(n) 

Vmi 

ysi 

ysi 

yi 

(s) 

TJ 

*(s) 

Z(s) 

•(n) 

w • 

_mi_ 

si 

■ 

SI 

Ci  • 

1 __ 

•T 

+ SS 

‘xfnf 

1 

+ S 

"•  (c)“ 
X1 

+ S 

uf 

■ 

,<n) 

• (c) 
yl 

Vf 

z!n) 

1 

•(c) 

_Z1 

wf 

td 

(4-1) 


This  expression  can  be  obtained  either  by  taking  the  expression  for  the  position 
of  a particle,  relative  to  the  origin  of  the  normal  earth-fixed  axes  and  resolved 
along  the  const ant- velocity  axes,  differentiating,  and  then  premultiplying  the 
result  by  S^S  to  refer  the  velocity  to  the  strip-fixed  axes;  or  by  adding 
together  the  various  relative  velocities  making  use  of  the  facts  that 
(of  Ref  4,  Appendix  A): 


(i)  If  some  axes  °axayaza  have  angular  velocity  <p^  r^  > relative 

to  axes  O^x^y^^  and  resolved  along  °axayaza  » then  the  linear  velocity  of  a 
point  relative  to  0.  resolved  along  the  a-axes  is 


where  ^ etc  are  the  standard  Euler  rotations  which  transform  to  the 
a-directions  from  the  b-directions , and  x^j^  etc  are  the  coordinates  of  0 
relative  to  0^  and  resolved  along  the  b-axes. 

(ii)  In  the  situation  described  in  (i) 


r (aTl 

• 

pb 

% 

9ab 

*ab 

(a) 

• 

qb 

eab 

(a) 

• 

l!b  J 

!ab 

(of  equation  (3-14)),  and 


(4-3) 


S, 

^ab  ^ab 

- S ST 
^ab  *ab 


(4-4) 


and 


(iii)  If  jx  y z J is  any  vector  and 


S any  axes  transformation  matrix  then* 


(4-5) 


In  equation  (4-1)  jx^  y^  is  given  by  equation  (3-8),  f being 

put  equal  to  zero  for  a particle  not  on  the  flap  portion  of  the  strip,  |<ji.  9.  <jj.| 
by  equation  (3-3),  and  |x.(n)  y<n)  z^n)J  by  equation  (3-2).  Substituting  theJe  1 
expressions,  and  making  use  of  (3-15),  (3-17)  and  (4-5),  we  have 


* (4-5)  is  of  course,  also  true  with  S replaced  1 ST  , ie 
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It  will  also  be  noticed  that  the  term  in  the  | } is  S^fR  (of  equation  (3-20)). 
The  square  of  the  resultant  speed  of  the  strip  reference  point  is  easily  seen  to 
be 


[uf  vf  wf]  uf  + 2 *r  - Ax.r  * + K 

if 

J :.(c)  • 


vf  1 +5tuf 


-A  <f>  + K q 

if  1 


4.2  The  forces  on  a strip 

4.2.1  Their  representation 

Consider  a vortex  whose  axis  lies  along  Oy  in  a frame  of  reference  Oxyz 
which  is  moving  with  linear  velocity  {u  v w},  referred  to  the  same  axes,  and  zero 
angular  velocity,  through  a fluid.  Then  the  force  on  that  vortex  is 


pA  0 
u 


p [ o o r fu 


o o o v 


- r o o w 


(4-9) 


referred  to  Oxyz  , where  T is  the  vortex  strength.  In  applying  strip  theory 
we  represent  the  strip  by  a set  of  bound  vortices*  whose  axes  are  parallel  to 


* Plus  also  a semi-infinite  layer  of  free  vortices  which  move  with  the  fluid  and 
so  have  no  force  exerted  upon  them. 
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°siysi  ’ ^he  stren8t^is  these  vortices  will  be  functions  of  the  boundary 

condition  and  so  of  the  velocity  of  the  strip  normal  to  its  surface.  We  will 

( s ) 

approximate  to  this  on  the  main  part  of  the  strip  by  w . at  points  on  0 .x  . ; 

mi  si  si 

and  on  the  flap  part  of  the  strip  by  the  velocity  normal  to  0g.y  ^ and  t0  the 
line  ( of  equations  (3-6)  and  (3-8))  given  by  y . - z . - 0 , ie  the  line 

81“  Sic 


r <sT 

“ 

X . 

si 

X,  . 

hsi 

y(s) 

■'si 

0 

Z(8) 

0 

_S1_ 

— _ 

+ (x  . - X.  .) 

sie  Tisi 


cos  6., 
if 


0 

- sin  6.- 
lf 


sin  6., 
if 


cos 


This  latter  velocity  is  therefore  the  last  element  of  P 


6i 


«(;) 

mi 

v(*> 

mi 

W^, 

Lmi  J 


(4-10) 

at  points  on 


(4-10),  which  to  the  desired  accuracy  is  {of  equations  (3-5)  and  (3-7)) 

[0  0 I] 


6if 


r.«i 

mi 

•f 

0 

0 

- 1 

ru(s)i 

mi 

(s) 
v . 
mi 

0 

0 

0 

v<8> 

mi 

w(?> 

mi 

1 

0 

0 

w<!> 

mi 

* [0  o 1]  P 


6if 


u<S> 

1 

v<8> 

l 

(s) 

w. 


* D 0 °1  p«.fVf 

if  if 


~ {\si  COS  6if  + (xsie  ' *hsi>}C°  1 °] 


(sT| 

i 

(s) 


(s) 
r . 

L i J 


~ - 

,T 

“ " 

uf 

f 

ql 

• 

vf 

• 

• 

_wf 

_qn 

- (x  . - x,  . )f 

sie  Tisi 


(4-11) 
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where  ju|S^  vf8^  w^|  and  |p£S^  r£8^}  are  8*ven  below.  It  will  be  seen 

( g \ , , # 

that  w . at  a point  on  0 .x  . on  the  main  part  of  the  strip  is  the  same 
mi  ^ si  si 

if 


mi  ” r~  "■  si  si 

expression  with  f and  6^  both  zero.  The  abbreviations  used  are: 


ps)l 

1 

r 

uf 

+ A S l Q F 

ql 

+ K 

• 

ql 

+ 

~-(c) 

X1 

+ A 

Uf 

♦ 

- A 

Xif 

‘il 

(s) 

# 

• 

* (c) 

0 

• 

v> 

i 

i 

vf 

• 

• 

• 

*1 

e 

(s) 

• 

• (c) 

• 

■ 

w: 

i 

wf 

_qn 

qn 

z] 

♦ 

4> 

(s) 

i 

(s) 


r<s> 

i 


qa  f 

♦ if 


m 

• 

ql 

+ S,fc 

♦if 

♦ 

• 

0 

• 

• 

• 

_qn 

♦ 

(4-12) 

(4-13) 


The  boundary  condition  is  therefore  made  up  of  certain  multiples  of  four  basic 
conditions  (at  the  points  described  above) : 


normal  velocity  of  air  = 1 


Xsie 


H(xhsi  ' Xsie) 

(xhsi  - Xsie)H(xhsi  ‘ xsie> 


where  H is  the  Heaviside  step  function. 
The  multiples  are: 


w.<S> 

1 


.(s) 


u<s)  sin  6.,  + wfS^(cos  - 1)  + qfS^x,  .(1  - cos  6.f) 
l if  l if  ’l  nsi  it 


and 


. , (us)  f (us) 

+ ,uif  cofl  6if  - wif  8 


in  6if)(«i  - 


6if> 


6i 


(us  denotes  unperturbed-strip-fixed  axes). 


I ■ 


L.  i'A l'  ” 


TM  Str  933 


TM  Str  933 


r 


- 

M 


i 


n 


19 


The  strength  of  the  vortex  will  therefore  be  a function,  for  a given  datum 
state,  of  these  four  scalars*  and  of  its  position  relative  to  the  wake  and  to  the 
points  where  the  boundary  condition  is  applied.  Thus  we  can  say  equivalently 
that  the  vortex  strength  is  a function  of 


u<S> 

1 


(s) 

w. 

1 


end  6. 

l l 


Thus,  for  small  perturbations  the  strength  of  a vortex  has  the  form 


r * rf  + ru(»<8>  - „<“*>)  * rv(v<8)  - „<*>)  . rq,{s>  * r{<5.  - «.f) 


(4-14) 


where  T^,  T^,  etc  are  functions**  of  u£fS^»  w£fS^  » an<*  ^if  • The  velocity 


of  the  vortex  relative  to  the  fluid  will  be  in  the  direction  of  the  boundary 
line.  The  enforcement  of  the  boundary  condition  ensures  that  there  is  no 
relative  velocity  normal  to  this  line.  The  force  on  the  vortex  is  therefore 
normal  to  this  line.  The  relative  velocity  of  the  vortex,  and  consequently  the 
force  on  the  vortex  will  be  a function  of  the  same  four  scalars  as  the  vortex 
strength  ( of  equation  (4-14)).  Thus  the  local  aerodynamic  force  at  a point  on 
the  strip  will  have  the  form 


C‘81 

1 

= Xif 

sin  5 . , 
if 

+ 

f<s> 

1 

0 

g<3> 
_l  _ 

cos  6.^ 

y.  sin5.r  Y-  sinfi...  sin6.,  y . r sin  S . c + y . , cos  6 

Aiu  if  Aiw  if  Aiq  if  Ai6  if  Aif  lij 


y.  cos  6.,  y.  cos  6.,  y.  cos  5.,  y . . cos  6 . c - x cos  6., 

Aiu  if  Aiw  if  Aiq  if  Ai6  if  Aif  if 


(s)  usH 
u;  - u;,  1 

l if 


w<S>  - w<“S> 
1 if 


.(■) 


. - <5-*  I 

L.  l if  — 1 


(4-15) 


* We  are  taking  the  view  that  the  strength  at  an  instant  is  determined  by  the 
boundary  condition  (and  a finite  number  of  its  derivatives  with  respect  to 
time)  at  the  same  instant. 
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where  6^  (in  the  trig,  functions)  and  (in  the  perturbation  part) 

are  put  equal  to  zero  off  the  flap  portion,  and  the  coefficients  y.  , y.  . v.  . 

Aiu’  *iw’  *iq* 

Xi6  may  *-nvolve  the  differential  operator  D . In  addition,  there  will  be  a 

leading  edge  'suction'  force,  as  will  be  seen  if  one  imagines  the  infinitely 

thin  leading  edge,  of  our  aerofoil  model,  as  the  limit  of  a rounded  edge.  The 

Kutta-Joukowski  theorem  tells  us  that  for  a two-dimensional  aerofoil  moving  with 

constant  and  uniform  velocity  the  overall  force  on  the  aerofoil  is  normal  to  the 

direction  of  motion.  The  leading  edge  'suction'  force,  which  has  only  a 

component  in  the  0 .x  . direction,  can  be  written  as* 
si  si 


X(2)  ♦ 
ifO 


'x< 

1 


(?)  x(?}  i<?> 

ixO  izO  Ai0O 


y(2)" 

ifiO 


(s)  (s) 

Ui  " Uif 

w<s>  - *<•> 
l if 

(s) 


6.  - 6.. 
l if 


(4-16) 


and  so  it  can  easily  be  shown  that 


- (2) 

ixO 


(4-17) 


(us) 

*if 

(us) 

if 


* r 


Y.  COS  6.^. 
1U  if 


-r 


X-  sin  6 . - + 0(D) 
Aiu  if 


(4-18) 


* The  (2)  superscript  indicates  two-dimensional  values,  and  the  circumflex  supra- 
script  has  been  used  to  indicate  the  affinity  of  these  coefficients  with  the 
body-fixed  axes  coefficients  of  Ref  1. 
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- £ xif  CO,  Sif  * -4^y  £ xi£  >in  S.£ 

if 

' {Z  Xi«  * Z Xi«  C0S  *4  " Z Xi«  *i£  * 0<D) 


(4-20) 


where 


^ indicates  a summation  over  the  main  part  of  the  strip,  and  ^ 


indicates  a summation  over  the  flap  part. 

The  overall  forces  on  a strip  are  ( of  (4-16))  written  as*,  referred  to  the 
strip-fixed  axes  (see  Fig  1) 


Xi  - X<f  ♦ Ig*  xg>  ig>  xgr  "u{8)-u^8)“  (4-21) 

yJs)  0 0000  wf8)  - wf“8) 

i l if 


zW  zfp  zf^ 

if  ix  iz  10  16 


6.  - 6., 

■ i if 


and  so  we  find  that 


Z»  ^<1 

xif  + Yj  Xif  C°8  6if 

'ix  = (Z  Xiu  + Z Xiu  C°9  6if  ' + °(D) 


{Z  Xiw  + Z Xi«  C°S  6 + °(D) 

{Z  Xi«  + Z Xi«  COS  6if  ' Z Xif  Sin  6if} 


(4-22) 


(4-23) 


(4-24) 


if  sin  6if>  + 0(D)  (4-25) 


* The  subscript  i has  been  used  to  indicate  a typical  strip  whereas,  in  Ref  I, 
, e.g.  was  the  force  due  to  unit  displacement  in  the  ith  mode. 
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r(2) 
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(us) 

if  Z<2> 


Uif 


(us)  "if 


(4-26) 
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*(2) 
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(us)  (us) 
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(U“S>) 


(us)  Z 


zf2)  + 0(D) 


Uif 
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(4-27) 


(us) 

Wif  ~ (2) 


Sir  z;r  - -ra-  zi£)  * 0<D) 


(4-28) 


uif 


w(us) 

Wif 


Uif 


-7 — v Z^2)  + 0(D) 
(us)  i<5 

‘if 


(4-29) 


Some  further  consideration  is  given  to  this  aspect  in  section  4.22  where  a 
further  relationship  (4-50)  is  obtained  (see  also  equations  (4-54)  to  (4-58). 
Similarly  the  overall  moments  on  the  strip  are  written  as,  referred  to  the  strip- 
fixed  axes*, 
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0 
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m(2) 

Mif 

IX 
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_ 
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1 if 

(s)  (us) 
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i if 
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6.  - 6., 
l if 


(4-30) 


There  will  be  no  contribution  to  these  from  the  leading  edge  'suction',  and  so 
we  find  that 


m(2) 

Mif 


M<?> 

ix 


y^,XifXsie  ^if^’Hisi 


si(cos  5if  - O 


■ EXiuXsie  “ (y  Xi^*h8i(cos  6if  - » 


(4-31) 


(4-32) 


* In  accordance  with  our  practice  the  symbols  for  the  moments  should  also  have  a 
subscript  s to  show  that  the  moments  are  about  the  origin  of  the  strip-fixed 
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"iz  y^xiw*sie  f ^ xiw)*hsi^COB  'J if  " 1 ^ (4-33) 

"if  ‘ ' £*i<,Ve  - (£  *i,)*hsi(<=»s  «i£  - '>  <4-3«) 

"if  ‘ ' EXi«Ve  -(E  *i«)*hsi(coB  6if  ' » 

*(Z  xif)’W  “n  ‘it  ■ «-35) 

The  strip  hinge  moment  can  be  obtained  as  a particularisation  of  this  obtained 
by  imagining  the  reference  point  (of  the  strip)  to  be  at  the  hinge  and  perform- 
ing the  summations  for  the  flap  part  only.  Thus  it  is 


where 
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*iu^Xsie  ^si^ 
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*i6^xsie  *hsi^ 


(4-36) 


(4-37) 


4.2.2  An  equivalent  form 

It  is  of  interest  to  express  the  aerodynamic  force  on  a strip  in  terms  of 
coordinates  which  describe  the  position  and  orientation  of  the  strip  relative  to 
its  position  in  the  datum  motion.  The  position  of  the  strip  reference  point. 


L 
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relative  to  its  position  in  the  datum  motion,  resolved  along  the  constant- 
velocity  axes  is  (of  section  3 and  equation  (3-2)) 


“(c) 
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xif 
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xif 
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(4-38) 


The  axes  transformation  matrix,  for  the  change  in  orientation  going  from  the 
unperturbed  position  of  the  strip-fixed  axes  to  the  strip-fixed  axes,  is 


ui 


*i  ♦ *if 


s.  s; 

♦i  *if 


s.  Vl 

*if  ♦ ^if 


(4-39) 


and,  so,  making  use  of  equations  (3-15)  and  (4-5),  the  standard  Euler  rotations 
which  achieve  this  change  in  orientation  are 
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(4-40) 


It  follows  therefore  from  equations  (4-12)  and  (4-13)  that 


. Am  • 


TM  Str  933 


Now,  since  is  zero,  S<j,^  is  given  by  equation  (3-24),  and  so 

cos  (0^  + - y£)  - sin  (6if  + ( 

cos  «|> sin(0£f  + 0f  - y£)  cos  <j> cos  (0££ 
0 0 

- cos  A.,.  sin  (0.r  + ©r  “ 
if  if  f 

cos  (0if  + 0f  - y£) 


+ sin  4>if 


0 0 
1 - sin  (0  . f + ©f  - Yf) 


sin  (0£f  + 0f 
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Incidentally,  the  angles  d> , 0^  + @£  - Yf  are  one  of  the  pairs  of  angles 
{of  Ref  3,  section  6.2)  which  can  be  used  to  specify  the  incidence  of  the  strip 
in  the  datum  motion.  That  is 


+ ®f  ” Yf  is  'incidence  magnitude' 


'if 


is  the  'incidence-plane  angle' 


Thus  (4-43)  is  an  expression  purely  in  terms  of  incidence  angles,  while  (4-46) 
is  its  transcription  in  terms  of  velocity  components  and  one  incidence  angle. 
The  latter  could  have  been  written  purely  in  terms  of  velocities  at  the  expense 
of  the  inconvenience  at  a J . 

From  equations  (4-21)  and  (4-46),  the  overall  forces  on  a strip  can  then 
be  written  as 
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(Similar  expressions  can  be  obtained  for  the  moments  and  hinge  moment,  from 
equations  (4-30)  and  (4-36),  but  we  will  not  detail  them  here.) 
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To  compare  with  the  form  taken  by  some  other  authors  ( eg  Baldock5)  it  is 
necessary  to  take  the  quasi-steady  form  - all  the  terms  containing  the  operator 
D in  the  above  rectangular  matrices  are  omitted  and  we  make  use  of  the  relation- 
ships (4-27)  to  (4-29)  - in  the  case  where  the  incidence  plane  angle  <t>^  of  the 
strip  in  the  datum  motion  is  small.  The  force  on  the  strip  then  has  the  magnitude 
(to  a first  order  approximation). 


lf  / I u(u»)| 

'"if 
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This  of  course  what  one  gets  from  two-dimensional  potential  flow  theory.  For 
example  for  a two-dimensional  aerofoil,  in  an  incompressible  fluid,  moving  with 
speed  (Vf  + <SV)  at  any  incidence  a , the  lift  is  proportional  to 

+ 6V)2  sin  (a  - a^)  where  is  the  no-lift  incidence.  Consequently,  the 

normal  force  on  the  aerofoil  is* 


Z « Zf  + Z^(u  - uf)  + Z* (w  - wf) 


C j-V^  sin  aQ  - 2V^  sin  (u  - u^) 

- vf)} 


+ cos  (w 


(4-51) 


where  C is  a constant,  and  so,  since 


uf  = Vf 


w^  = 0 


(4-52) 


Z*u,  + Z*w- 
x f z f 


- 2CV^  sin  (*q 


2Z . 


(4-53) 


The  relationship  (4-50)  did  not  appear  from  the  reflections  of  section  4.2.1 
since  there  we  took  account  only  of  our  knowledge  of  the  direction  of  the 
resultant  force  on  the  aerofoil  rather  than  its  magnitude.  If  (4-50)  is  com- 
bined with  equations  (4-26)  to  (4-29)  we  have  the  relationships 
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* We  are  taking  = 0 for  simplicity,  but  this  involves  no  loss  in  generality. 


The  same  relationship  (4-53)  is  obtained  for  any  . 
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(of  equations  (3-14),  (3-24)  and  (4-7))  and  so,  from  (4-21),  the  overall  forces 
on  a strip,  referred  to  the  strip-fixed  axes,  are 


From  (4-36)  we  see  that  the  hinge  moment  B.  has  the  same  form  as  the  pitching 
(s)  1 

moment  - all  the  M coefficients  are  merely  replaced  by  B coefficients. 

The  final  forms  of  the  above  expressions  can  of  course  be  used  even  when 

(2)  (2)  (2) 

one  is  not  using  two-dimensional  strip  theory  - X.,  , Z.,  , M.,  (with  the  2 

Spy  It  it  It 

replaced  by  a 3),  and  the  row  vectors  X.  etc,  being  given  appropriate  three- 

iq 

dimensional  values. 

4.3  The  overall  forces  and  moments  on  the  aircraft 


The  overall  forces  on  a strip,  referred  to  the  constant-velocity  axes, 
are,  using  the  relationshios  (3-15)  and  (4-5)  along  with  equation  (4-60) 
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Similarly  the  overall  moments  on  a strip  about  the  strip  reference  point*, 
referred  to  the  constant-velocity  axes,  are 
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(4-63) 


* In  accordance  with  our  practice  the  symbols  for  the  moments  should  also  have  a 
subscript  a to  show  that  the  moments  are  about  the  origin  of  the  strip-fixed 
axes  but  this  has  been  omitted  to  avoid  confusion  with  the  subsequent  use  of  a 
subscript  s to  denote  structural. 
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4.4  The  generalised  forces 

The  generalised  aerodynamic  forces  could  be  determined  using  the  expressions 
derived  in  Ref  1 , and  generalised  in  Appendix  A,  which  are  in  terms  of  the 
overall  and  local  forces  referred  to  the  constant  velocity  axes  ( of  equations 
(A-18)  and  (A— 19)).  However,  for  the  particular  semi-rigid  and  aerodynamic  models 
being  considered  in  this  paper  it  is  more  convenient  to  proceed  rather  differently. 


The  linear  velocity  of  the  strip  reference  point,  and  the  angular  velocity 
of  the  strip-fixed  axes,  both  resolved  along  the  strip-fixed  axes,  are,  for  tha 
ith  strip,  respectively  (of  equations  (4-1),  (4-3),  (4-41)  and  (4-42)) 


1 

1 

(s) 


1 

• T 

+ SS 

(n) 

x . 

l 

+ S 

‘(O 

X1 

+ s 

uf 

y(n> 

} i 

•(c) 

yl 

vf 

1 

2<n) 

1 

• (c) 
Z1 

wf 

and,  by  a well  known  relationship  ( of  for  example,  Ref  4) 

= Q* 


(sf| 

i 

(s) 


r<S> 

l 


ui 


ui 


e . 

Ul 


Ip  . 

Ul 


(4-77) 


(4-78) 


where  the  rotations  $ ^ etc  satisfy  equation  (4-39),  Q<j>u^  is  defined  by 
equation  (3-14),  and  it  can  be  shown  that,  using  (4-39) 
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* This  expression  can  be  deduced  considering  the  angular  velocity  of  the  strip- 
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A (s)  = S*  . 

p>  ' YU1  U1 
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•T  *T  T 
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(4-79) 


constant 


The  position  of  the  strip  reference  point,  relative  to  the  origin  of  the 
mt-velocity  axes,  resolved  along  the  strip-fixed  axes,  is  given  by,  say 


and  so 
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(4-80) 


(4-81) 


Consequently,  from  the  work  done  in  time  <5t  , going  to  the  limit  6t  = 0 , we  find 
that  the  virtual  work  done  by  the  forces  on  the  strip  in  small  displacements 

6ql»  *•*»  5qn+6  1S 
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The  forces  and  moments  which  appear  in  the  strip  contribution  to  the  generalised 
force,  equation  (4-82),  are  given  by  equations  (4-60)  and  (4-61 ),  and  so  using 
these  equations,  along  with  (4-89)  and  (4-91),  we  have,  writing 
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(4-92) 


TM  Str  933 


TM  Str  933 


F 


47 


and 
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where  £ is  the  lower  triangular  matrix  formed  from  the  elements  of  {$  9 i|/|  viz 
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(4-95) 
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(2)  (2)  (2) 

In  addition  Y:r  , L.'  , N;  ' are  each  zero  and  so 
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(4-98) 


The  second  and  third  rows  of  submatrices  in  the  matrix  [Q^j]  are,  as  expected 
(of  Ref  1,  equations  (74)  to  (78)  and  Table  3),  closely  related  to  the 
coefficients  in  the  expressions  (see  equations  (4-64),  (4-72),  (4-74)  to  (4-76) 
and  (4-62))  for  the  overall  forces  and  moments.  Making  use  of  the  general 
relationships  (4-96),  (4-5),  and 
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we  find  that 
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- (see  Equation  (4-93)) 
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These  expressions  for  the  aerodynamic  contributions  to  the  generalised  forces 
(equations  (4-93),  (4-94)  and  (4-100))  are  not  restricted  to  the  two-dimensional 
strip  theory  aerodynamic  approximation.  As  remarked  in  section  4.23  they  can  be 
used  in  general  provided  the  appropriate  meaning  is  given  to  the  various  aero- 
dynamic coefficients. 


4.5  Strip  interference 

An  important  factor  affecting  the  dynamical  behaviour  of  an  aircraft, 
particularly  when  it  is  moving  largely  as  a rigid  body,  is  the  aerodynamic  inter- 
ference between  wing  and  tail.  This  is  a particular  manifestation  of  the  fact 
that  there  is  in  general  such  interference  between  any  two  parts  of  the  aircraft, 
ie  in  our  representation,  between  any  two  strips.  So  far,  in  our  analysis,  we 
have  largely  ignored  such  interference.  The  various  formulae  that  we  have 
obtained  can,  as  we  have  pointed  out,  be  used  with  three-dimensional  values  for 
the  aerodynamic  coefficients,  and  so  full  account*  can  be  taken  of  the  inter- 
ference effects.  For  example,  the  expression  (4-21)  for  the  overall  forces  on  a 
strip  referred  to  the  strip-fixed  axes  could  be  extended  to  include  terms 


* Provided  of  course,  that  one  has  a useable  aerodynamic  theory  for  the  whole 
configuration.  Refs  8 and  9,  for  example,  go  a good  way  towards  providing  this 
for  the  subsonic  case. 


proportional  to  the  perturbations  of  the  velocities  and  control  angle  of  the  other 
strips  - one  or  more  of  them  - and  the  perturbations  of  their  relative  positions 
and  orientation.  Thus  we  can  write 
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(4-101) 


and  similar  expressions  for  the  moments  and  hinge  moment  on  the  ith  strip.  In 
these  expressions  Ax^  etc  are  measures  of  the  change  in  the  relative  position 
and  orientation  of  the  ith  and  jth  strips.  The  coefficients  xf?}  etc  will  be 

1XJ 

functions  of  the  position  and  orientation  of  the  ith  strip  (ith  strip-fixed  axes) 

relative  to  the  other  strips  (in  the  other  strip-fixed  axes)  in  the  unperturbed 
• • • A (s ) 

state.  The  coefficients  in  the  last  term,  X. . etc,  can  be  taken  as  constants  - 

...  . .(s  1JX  . . 

the  coefficients  in  the  middle  term,  X^(  etc,  will  in  general  be  differential 

operators  - since  any  differential  operator  terms  can  be  included  in  the  other 

coefficients . 

4.5.1  Using  two-dimensional  theory 

Equation  (4-101)  indicates  the  possibility  of  making  use  of  two-dimensional 
theory  for  interfering  surfaces  when  one  hasn't  an  adequate  three-dimensional 
theory.  It  is,  however,  difficult  to  do  this  in  general.  A rather  crude 
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approximate  way  of  doing  this  is  therefore  suggested  below.  In  certain  circum- 
stances it  will  become  a good  approximation,  and  in  any  case  it  should  be  in  the 
right  street. 

If  the  y strip-fixed  axes,  of  the  ith  and  jth  strips,  were  parallel,  and 
the  x strip-fixed  axes  were  coplanar,  in  the  datum  state,  and  the  datum  flight 
path  direction  was  parallel  to  that  plane,  then  we  could  consider  these  strips 
as  part  of  a two-dimensional  configuration  and  obtain  the  aerodynamic  interference 
accordingly.  However,  this  is  an  exceptional  situation.  The  way  the  aircraft  has 
been  divided  into  strips  ensures  that  the  last  two  conditions  are  satisfied  when 

yif  * yjf  ■ «-,02> 

but  to  satisfy  the  other  condition,  we  must  have 
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(4-104) 

It  is  therefore  suggested,  as  a rough  approximation,  that  one  should,  for  any 
pair  of  strips  which  satisfy  (4-102),  obtain  (two-dimensional)  interference 
aerodynamic  forces,  assuming  (4-103)  was  true,  and  then  multiply  these*  by 
cos  (d> ^ - 4k)  . Thus  (4-101)  is  rewritten  as  (c/  equation  (3-3)) 


* The  difference,  between  a force  in  the  isolated  and  tandem  configurations,  is 
multiplied  by  cos  (4>^  - $ . ) and  then  added  to  the  force  for  the  isolated 
configuration.  J 
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(4-105) 


where  the  jth  strip  is  one  which  satisfies  equation  (4-102)  - we  are  assuming, 
as  will  almost  certainly  be  the  case,  that  there  is  only  one  such  strip  - the 
subscript  (i  - j)  denotes  the  difference  between  the  values  of  the  modal  function 
at  the  ith  and  jth  strips,  and 


= (xjn)  - ,<»>)  cos  0.  - (.<»>  - z<n))  sin  0.  J 


- (*jf  - x.f)  cos  6.f  - (Zjf  - z.f)  sin  e.f 


(4-106) 
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The  use  of  equations  (3-2),  (3-3),  (3-7)  and  (4-59)  will  enable  one  to  write 

(4-105)  entirely  in  terms  of  the  generalised  coordinates.  The  expressions  for 

the  moments  and  hinge  moment  at  the  ith  strip  will  be  similar  to  equation  (4-105). 

(2)  A (2)  *■  (2) 

The  interference  coefficients  in  the  expressions,  such  as  X..,,  X...,  Z...  etc, 

r ljf  1X1  1XJ 

will  be  functions  of  x^ J , z^J , 0*J  , u^,  w£f>  ^if»  ^jf*  ^n^er  a^a  (of  Fig  2), 
where 
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(4-111) 

To  avoid  ambiguity,  the  coefficients  with  subscripts  ijx,  ijz  or  i j 0 


A (2 ) 

(eg  X. . ) are  taken  to  be  constants. 
a ijx' 

operators . 


The  other  coefficients  may  be  differential 


For  some  the  significance,  in  some  respects,  of  these  descriptions,  such 
as  (4-105),  for  the  forces  and  moments  on  a strip,  may  be  illuminated  if  we 
consider  the  'purely  two-dimensional'  case,  when* 
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* From  our  specification  of  the  datum  motion  (section  2)  and  the  division  into 
strips  (section  3),  i|k  ^ and  v^  are  always  zero. 


and  similarly  at  the  jth  strip,  and  the  perturbations  are  expressed  as  displace- 
ments relative  to  an  inertial  frame.  No  loss  in  generality  is  achieved  by  taking 

this  frame  to  be  the  unperturbed  strip-fixed  axes  (for  each  strip)  and  the  body 

(c) 

freedom  displacements  Xj  , <j>,  to  be  zero.  Then  (of  equations 

(4-59),  (3-2),  (3-3),  (3-7),  (3-15)  and  footnote) 
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Also 
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and  so,  for  example,  the  forces  on  the  ith  strip,  in  the  direction  of  the 
unperturbed  strip-fixed  axes,  are: 


TM  Str  933 


I 


55 


.r-x!2>.x<?>D  -i'2).x<;;» 

ljx  ixi  Xjz  iti 


zf?*  -Z<?>  ♦ zfjjp 

ijx  1x1  Ijx  iil 


Mil - (MU * MJMil - *1 JxtU) * Mils "■  #if * uf *ineif>  * Mil s “•  8if -f 9i£>  * MSI0  K,i 

0 o 

- 'Iff  - (Mil  * - -I1*"!)  - — 9if  ♦ > »( ■ *1, > ♦ Mils  »•  *1,  ■-•if)  • Mil0  M.i 


(us)  (us) 
1 '*if 

(us)  (us) 
li  ’“if 


(Mil  •••  MJ  - Mil  •*•  Mj)  ♦ Mil0  (Mil  «•  Mj  * Mil  MJ)  • Mil0 

0 0 

(Mil — M1  - Mil  “■  Mj)  * M‘1]0  (Mil  •••  Mj  * Mil  •*•  MJ)  * Ml]° 

Mi#  * M«i<vf c°* *jf  *“f ,in *jf > • Mil s •••  'if - *f ,io 'if*  * Mji°  MS] 

0 0 

Ml«-Mii(“fco"jf*"f,ia'if  ♦M!i<“f,»,if-f'‘«v*MiiD  Ml] 


(us)  (us) 

i 'i« 


This  can  be  written  as 

r„(us)i  ^ r.(2) 


x<us) 

1 


zfus) 

1 


(4-115) 


x:  ' x:  ' 

ixk  lzk 


k=i,  j 


(2)  (2)  (us)  _ (us) 

i9k  16k  \ 

o o 


7 V'-  / 7'‘-' 

i6k  i6k 


0.  - 
k kf 


<5,  * 

k kf 


in  the  ’purely  two-dimensional  case’,  where,  with  j 4 i 


(2) 

« (2) 

= - X.  . + 

- (2) 
x):.d 

ixi 

IJX 

1X1 

(2) 

- (2) 

= X.  . cos 

ejj  - 

xfP  sin  + xf^ 

1JZ  f 1XJ 

jxj 

Ijx 

f 

(4-116) 


(4-117) 


(4-118) 


56 


(2)  . 

. -x<2> 

+ X^??D 

(4-119) 

‘izi 

ljZ 

1Z1 

(2)  . 

= cos  0.^  + 

y(2)  . 

X.  . sin 

6fj 

a(2) 

+ x-*-D 

(4-120) 

izj 

ljZ 

t 

1JX 

f 

i*J 

.(2) 

‘i0i 

- Z<2>- 

ijf 

(*$  ♦ 

z^<2)  - 
f ijx 

i: 
■ xf 

ixf2>) 

ijz  / 

- 

(w^  cos 

0if  + Uf 

sin 

0.f)xf2^  + 

if  1X1 

(u,.  cos  0.r  - wc  sin  0 . , ) X . • . 
f if  f if  izi 

* (2) 

+ X.f.D 
101 


(4-121) 


y(2) 

Xi9j 


y(2) 

xije 


* (2)  . « (2) 

(wf  cos  6jf  + uf  sin  0jf)X££j  + (uf  cos  0jf  ~ wf  sln  9jf^Xizj 

(4-122) 


+ X^D 
i0J 


X 


(2) 

i6i 

c<2> 

i<5j 


y(2) 

Xi6i 


= X 


(2) 

i6j 


(4-123) 

(4-124) 


and  similarly  for  the  Z coefficients  with  in  particular 

z<2>  ■ - x!2>  - (z <2>  * 4iz<2>  - 

101  ljf  \ 1J0  f ljx  f ljz/ 


^(2)  *(2) 
- (w f cos  9.f  + uf  sin  0if)Z\(  + (uf  cos  0.f  - wf  sin  0if)ZiH 


+ 

i0i 


(4-125) 


These  equations  enable  one  to  derive  the  ' circumflexed'  coefficients  from  the 
' uncircumf lexed'  coefficients.  If  we  introduce  the  notation* 


e(  ) = steady*  part  of  ( ) 

3(  ) = (1  - e)(  )D_1 


(4-126) 

(4-127) 


* Thus  for 


F = Fq  + F,D  + F2D 


e(F)  = F 


0 


3(F)  - Ft  + F2D  . 


I 
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and  write,  j j 1 


X*  . 
101 


X*. 

10J 


z* 

i0i 


Z* 

iej 


(2) 

X;Q.  + (w,  cos  0 . , + u,sin0.r)3 
i0i  f if  f if 


(Xixi)  ' (ufc°s0if  - wf  sineif)8(Xizi) 

- z iUxW)  + xjjE/x<2>)  - Z<2> 

f \ lxi/  f y iz i/  ijf 

Xi0j  + (wf  COS0jf  + uf  sin0jf)3(Xixj)  " (uf  cos0jf  - vf  sin0.f)a(x£)) 


(4-128) 


(4-129) 


Z<2>  ♦ 
101 


(wfcos0.f  + uf  sin0if)9(z^)  - (Ufcos0.f  - wf  sin  0.f ) a(z<2>) 

- z^e(zf2})  + x^e(zf2^  + X^£ 

f \ ixi/  f y izi/  ijf 

if)3(zi«j)  ' <uf  cos8j£  • ”£  sin  8j£)8(ZS) 


(4-130) 


(2) 

Z'  . + (w  cosfl,,  + u,sin0., 
10J  f jf  f jf 


(4-131) 


Then  we  have 


~x<2) 

ljX 

x<2) 

ijz 

$(2) 

Xije 

= e 

-x<2> 

1X1 

-x<2> 

1Z1 

-**ei 

0 

0 

0 

0 

0 

0 

z.(2) 

ljX 

z<2) 

1JZ 

4(2) 

ijej 

-z<2> 

1X1 

-z<2> 

1Z1 

-z* 

Ziei 

_ 

xf2J  cos  0*j  + xf2?  sin  0^j  X^2}  cos  0^  - X^2}  sin  0^ 

ixj  f izj  f izj  f ixj  f 


ZS  =°s  «£j  * 41]  si"  e£j  ZUj  *°s  9fj  - ZS  8i”  8£j 


X* 

iej 


and 


1(2) 

ixk 

0(2) 

iSk 

i<l> 

10k 

= 3 

" (2) 
ixk 

Y(2) 

izk 

X* 

i0k 

0 

0 

0 

0 

0 

0 

4(2) 

ixk 

4(2) 

izk 

z<?> 

i0k 

_4ll 

Z-2fc 

izk 

z* 

i6k 

(4-132) 

(4-133) 


These  two  equations,  along  with  (4-123),  (4-124)  and  their  partners,  give  all  the 
'circumf lexed'  coefficients  required  in  (4-105).  Similar  equations  will  give  the 
coefficients  appearing  in  the  expressions  for  the  moments  and  hinge  moment. 

For  this  'purely  two-dimensional'  case  the  moment  about  the  strip  reference 
point  referred  to  the  unperturbed  strip-fixed  axes  is 


and  so  the  moment  about  the  origin  of  the  unperturbed  strip-fixed  axes  is 


(4-134) 
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(s) 

while  writing  M)  ' in  the  form  of  (4-105)  and  using  (4-113)  and  (4-114)  gives, 


after  some  analysis,  the  relationships: 


i 


(2) 

ijx 


= - eL*  . M*  . M?fll 

ljz  1J0  1X1  1Z1  10lJ 

= e ImP?  cos  0^  + sin  0^  lif2?  cos  0^  - M^2}  sin  0^  M*  1 
1XJ  f 1ZJ  f 1ZJ  f 1XJ  f 10J 


[SUk  “ilk  “Ilk]  • 8[“Lk  "izk  “iekj 


where 


and 


fi(2>  _ M(2) 

Mi6k  " Mi6k 


* 

Mixk 


Mizk 


1X1  ljf 


M(2) 

Mixk 


mP?  - 

1Z1  ljf 


= M 


(2) 

izk 


(4-137) 

(4-138) 

(4-139) 

(k  - i) 

(k  * i)  (4-140) 

(k  “ i) 

(k  + i)  (4-141) 


M*  , 
10k 


- m<2>  ♦ 
101 

(wf 

cos 

+ 

M-t 

•r-l 

CD 

uf 

sin 

eif>8(MS) 

- 

(uf 

cos 

9if  “ 

wf 

sin 

eif>s(Miii) 

- 

[*$ 

+ e^l 

<<21 

IX] 

!)]- 

• 4$  - •( 

>sj 

(k 

« i) 

- mP  + 

i0k 

(Wf 

cos 

9kf  + 

uf 

sin 

v9(MS) 

- 

(Uf 

cos 

9kf  " 

wf 

sin 

ekf>8(“llk)  • 

(k  * i) 

(4- 

•142) 
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E; 


The  corresponding  relationships  for  the  hinge  moment  coefficients  are  easily  seen 


to  be 

|sf2> 

1 1JX 

s<2> 

1JZ 

fi(2) 

Bij0 

= - cIb^  B f 2 ^ 

1X1  1Z1  101J 

= e Jb^2^  cos  0^  + B^2?  sin  0^ 

1 1XJ  f 1ZJ  f 

B^2?  cos  0^ 
izj  f 

-Bf2>sine^  B*  . 

1XJ  f I0jl 

(4-143) 

fs 

;(2) 

izk 

B.(?) 

i0k 

■ a[Bi*t  BS  Bk] 

(4-144) 

S(2) 

Bi6k 

■ B«k 

(4-145) 

where 

B* 

i0k 

(2) 

= Bj_gi  + (wf  cos  0^  f + u£  sin 

9if>S(Bl*l) 

- (uf  cos  0^  - wf  sin 

9if>s(B"i) 

- «?« 


(b^2})  + x^e^2}) 

\ 1x1/  f \ izi/ 


(k  = i) 


= Bi0k  + (wf  cos  9kf  + uf  8in  9kf>3  (Bixk) 

- (uf  cos  ek£  - wf  sin  O^fc2^  (k  i i)  (4-146) 

^•5.2  And  in  terms  of  the  generalised  coordinates 


From 

equations  (3-2) 

and 

(3-3),  it 

Ax 

cos  9., 
if 

0 

- sin  9 . , 
if 

Az 

< 

sin  0 . c 
if 

0 

cos  0., 
if 

A0 

0 

0 

0 

(K)._i  ♦ 


0 0 
0 0 
I 0 

ij 


-x1J 

xf 


(F)j 


(F>i 


Lqnj 


(4-147) 
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Equation  (4-59)  gives  the  other  coordinates  that  we  require  in  terms  of  the 
generalised  coordinates.  The  coefficients  in  equation  (4-105)  are  given  by 
equations  (4-123),  (4-124),  (4-132)  and  (4-133).  Putting  all  these  together  we 
obtain  the  following  expression  for  the  overall  forces  on  the  ith  strips 


J 

!i 
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The  interference  coefficients  ( ie  all  the  ones  with  a triple  subscript)  in  the 
above  expression  will  be  functions  of  x^ , z^J , 0^,  u£f*  wjf»  and 

(cf  equations  (4-109)  to  (4-111).  Similar  expressions  can  be  obtained  for  the 
moments  and  hinge  moment  on  a strip,  to  take  some  account  of  the  aerodynamic 
interference  between  strips,  using  the  relationships  (4-137)  to  (4-139)  and 
(4-143)  to  (4-145). 

5 THE  OTHER  CONTRIBUTIONS  TO  THE  GENERALISED  FORCES 

The  contributions  to  the  generalised  forces*  from  the  other  forces* 
(gravitational,  structural,  etc)  acting  on  the  aircraft  can  be  obtained,  as  is 
most  convenient,  either  by  obtaining  expressions  for  the  translational  forces, 
moments  and  hinge  moment  on  a strip  and  proceeding  as  we  have  done  for  the  aero- 
dynamic forces  (e/  section  4.4  and  in  particular,  equations  (4-82),  (4-89)  and 
(4-91)),  or  by  substituting  the  expressions  for  the  local  forces  and  modal 
functions  in  the  formula  obtained  in  Ref  1 (section  6.1)  and  extended,  for  our 
form  of  deformation,  in  Appendix  A of  this  paper  (equations  (A-12)  and  (A-13). 
With  the  latter  approach,  the  modal  functions  required  are  (see  equations  (3-20) 
and  (A-7)): 


R = K - S*  A . . Q , F + P* 

♦if  x(u^  ♦if  *if  6if 
sif 


0 0 1 

x . - x,  . 

sie  nsi 

0 0 0 

y . 

sie 

- 1 0 0 

z . 

sie 

(5-1) 


at  the  ith  strip  , 


M„ 


SI  \ F 

♦if  *if 


(5-2) 


'if 


r 

•« 

-A  , .Q.  F + 2p"F 

x(YS>  ^if  6if 

< Slf 

0 

0 

1 

Xsie  ^Sisi 

fT 

0 

0 

0 

^sie 

- 1 

0 

0 

z . 

sie 

_ 

(5-3) 


* The  word  force,  without  the  adjective  effective,  is  to  be  understood  as  meaning 
applied  force.  The  distinction  is  between  applied  forces  and  reversed  effective 
forces  (of  section  6). 
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ST  PT 

♦if  6if 


(5-4) 


T T 

- s!  p;  a 


6.  "x  . -x,  . 
if  if  sie  hsi 


(5-5) 


U = cos  0ifS^ 


ifAx(uS>(Q*if) 


(5-6) 


With  the  former  approach,  that  is  proceeding  as  in  section  4.4,  if  typical 
forces,  moments  and  hinge  moment  on  a strip  are,  referred  to  the  strip-fixed  axes, 


IITIJ  _ 

z 


x(us)  + 

if 

y(us> 

if 

j(us) 

if 


ik  qk 


(5-7) 


PTO 

-Z  E» 

5<“8>  k"  5.. 

if  ik 

5^s)  5.. 

if  ik 


ik  ^k 


(5-8) 


B.  = B.,  + 

l if 


n-fo 


(5-9) 


then  the  contributions  of  these  forces  to  the  generalised  forces  are  found  to  be 


n+t> 

* Z 5jk\ 


(5-10) 
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referred  to  the  constant-velocity  axes  where  is  the  last  column  of  S^, 

f f 

and  the  vector  {<!>  4^}  is  given  by  equation  (2-2).  Thus 


sin  4*^  cos 


cos  4>  cos  ©j 


cos  0, 


(5-14) 


Consequently  the  local  gravitational  force  referred  to  the  strip-fixed  axes  is 


(5-15) 


and  the  overall  forces,  moments,  and  hinge  moment  on  a strip  are,  respectively, 
making  use  of  equations  (3-15)  and  (3-3),  and  defining  iik  as  the  mass  of  the 
ith  strip  ( of  equation  (6-17)). 


mi8Vsl*e 

1 f 


mig  S( 


’u\  + [s»i 


A S':  Q.  F 0 S 
if  \ 


, A8.  ql 

'if  \ 

I • 


(5-16) 


L, 


g I 4mv J V 

\atrip  ai j 


S i**  [a  s 
‘iO  *if  *f  [ aiO  *i 


A,  SI  <h  F - A P.  S.  £*  fT  0 A 
if  **f  *il  *lf  ‘ if  ai3  4if  *if  *f 


■A,\]h ' 

l5n+&_ 


(5-17) 


B^>  - g[0  , 01  J k Si 

Vflap  xsi  \si 


n.g  - aT  P S l + f-  aT 

1 < 13  4u  *i  l 1 


P S A sl  Q p - aT,P.  s i*  fT  0 

3 4if  *if  l<b  11  5<*  *;r 


aT  P S,  A,  irq, 
13  4if  ♦«  ' 


(5-18) 


where,  when  the  flap  angle  is  zero,  referred  to  the  strip-fixed  axes 

the  strip  eg  is  at  |xio  ® oj 
and  the  flap  eg*  is  at  j^si  + xi2  0 * 


and  where  m. , iniik  are  respectively  the  mass  of  the  strip  and  of  the  flap  part 
of  the  strip,  and 


* For  an  unmassbalanced  flap  one  will  expect  x^  t0  be  negative, 
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ail 


= ,n'.x.„ 

1 i2 


i3 


TT.X.„ 
i i2 


0 

0 

1 


(5-19) 


(5-20) 


3i0 


- TT 


• X 
1 


i2 


0 

0 


il 


(5-21) 


In  this  case,  it  is  clearly  much  simpler  to  use  the  formulae  of  Appendix  A 
(equations  (A— 1 2 ) and  (A- 13))  rather  than  carry  out  a similar  analysis  to  that  of 
section  4.4.  Using  (5-13)  we  then  find  that  the  gravitational  contribution  to  the 
generalised  force  in  the  jth  degree  of  freedom  is 


n+6 

“ Gj  = “ Gjf  " Yj  GJkClk  (5~22) 

k=l 


— 

r — v / 

where 

G,f 

= “ 8 

) m.  KT  + FV  AS-  faT  P S ) 

4"  l\  *i£  ai0  *if  13  Sif  *u) 

• 

ml 

Gn+6 , f 

Y m.A 

X x.f0 

with 


Xif0 

= 

_Xif" 

yif0 

yif 

Zif0 

Zif 

a 

if 


iO 


(5-24) 


and  I a 3 x 3 unit  matrix,  and 
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H* 


5.2  The  structural  contribution 


In  this  case  it  is  simplest  to  use  equations  (4-91),  (4-89)  and  (4-82) 
where  the  forces  and  moments  in  the  latter  equation  are  the  structural  forces, 
moments  and  hinge  moment  on  a strip.  These  forces  will  be  independent  of  the 
rigid  body  perturbation  of  the  aircraft  and  so  we  can  write  them  respectively  as 


(5-31) 


an + E 


Lsik  qk 


(5-32) 


»<!>  + t B • 

sif  [_j  si 


sikqk 


(5-33) 


Now  the  structure  cannot  exert  any  overall  force  or  moment  on  itself  and  so  it  is 


easily  shown  that 


I « 


(5-34) 


j 
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and  the  submatrix  E is 

qq 


E = \ 

qq  Lj 

i 

kV 

♦if 

i 

r 

X ...  X .“ 

sil  sin 

Y 

“A  , vQ,  FI 

si  1 

Z 

- 

sil 

m 

T 

+ F 

< 

r t 

L • • • L • 

sil  sm 

M 

• 

- cos  0 . c t F 

K,‘S) 

sil 

N 

a 

_ Sil 

\ 


+ f IV?}  ...  bV]  . (5-41 

I si I sin) 

This  all  appears  very  simple,  but  of  course  quite  a lot  of  work  may  be  involved 

(s) 

in  getting  the  structural  forces,  Xg^  etc,  on  a strip.  If  the  strip  can  be 
considered  as  a slice  of  a beam  then  these  forces  and  moments  are  the  increments 
in  the  shearing  forces,  bending  moments,  etc,  going  from  one  side  to  the  other 
of  the  strip. 


5 . 3 The  propulsive  contribution 


In  Refs  1 and  2 a very  simple  model  of  the  propulsive  force  was  used. 
Improvement  of  that  model  is  being  considered,  amongst  other  things,  in  another 
paper^  currently  being  written.  We  will  therefore  in  the  present  paper  stick 
to  the  simple  model  though  for  convenience  it  is  not  quite  the  same  model  as  that 
used  previously ^ . Thus  we  assume  that  the  propulsive  force  acting  on  a strip 
is  such  that  it  has  constant  components  in  the  direction  of  the  strip-fixed 
axes,  and  similarly  for  the  propulsive  moment  about  the  strip  reference  point. 
That  is,  we  assume 


X(?T 

P1 

= 

"(sf 

XPif 

Y(s) 

Y(S> 

pi 

pif 

z(s) 

z (3l 

pi 

pif 

(5-42) 


We  also  assume  that  there  is  no  propulsive  force  of  moment  on  the  flap  part  of 
any  strip. 

The  propulsive  contribution  to  the  generalised  force  in  the  jth  degree  of 
freedom  is  then  easily  seen  to  be 


n+6 


k=l 


(5-44) 


where  {.of  section  4.1) 


(5-45) 


(5-46) 
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and  the  overall  propulsive  forces  and  moments,  about  the  origin  of  the  constant- 
velocity  axes,  and  referred  to  those  axes  are  (of  section  4.3) 


74 
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(5-52) 
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With  our  modal  matrix  R , equation  (3-20),  it  is  easily  shown  that 


P 

q 


(5-55) 


where  the  subscript  0 indicates  that  the  expression  is  evaluated  for  the  strip 
containing  the  reference  point.  If  therefore,  the  only  propulsive  forces  acting 
on  the  aircraft  are  acting  on  the  main  part  of  this  strip,  the  submatrices  in 
(5-46)  have  the  particular  values: 


KTA  P - cos  0 . f FT£ 

V q lf  / 


Q*'1  S. 

*if  *if| 


pf 

V 

LNPfJ 


- A 


Pf 

Ypf 

LZP£J 


'll 


(5-56) 


P 

<t><7 


" f(K)o 

Pf 


(5-57) 


(5-58) 


(5-59) 


(5-60) 


Thus,  in  this  particular  case,  when  the  present  assumption  and  that  of  Ref  1,  as 

regards  the  propulsive  forces,  become  identical,  all  these,  except  for  P 

qq 

are,  as  expected,  the  same  as  the  expressions*  obtained  in  Ref  1,  Table  3.  The 
difference  in  the  P , as  in  the  other  contributions  to  the  generalised  forces, 
arises  from  the  fact  that  the  expressions  for  the  deformations,  here  and  in  Ref  1, 
only  agree  to  first  order  in  the  generalised  coordinates  ( of  Appendix  A). 


* On  the  main  part  of  the  strip 

R = 


containing  the  reference  point 


The  velocity  of  the  constant-velocity  axes,  resolved  along  the  strip-fixed 
axes  for  the  ith  strip,  is,  using  (3-15), 


Consequently,  the  velocity  of  a particle  relative  to  the  constant  velocity  axes, 
and  resolved  along  the  strip-fixed  axes,  is,  from  (4-6), 
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The  kinetic  energy  W , of  the  system,  relative  to  the  constant  velocity  axes, 
is  given  by 


(6-3 


and  so  it  is  easily  seen  that  the  generalised  effective  forces  are 
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6mA 


(us) 


strip  flap 

In  passing,  we  note  that  this  means  that 

I 


strip  sif 

sie  ’Sisi' 


Z6mA , , 


m.  A 

1 aio 


m.  A 
1 ai1 


6mA  = - I.n  ~ m.  A 

x • 10  1 , *hsi 


strip 


sie 


A P.  - A 
6 . , a . . 6 . - a. 
if  ll  if  ii. 


(6-10) 


(6-11) 


+ (Pt  A P.  - A >A 
1 6if  ail  6if  ill  ^isi 


* (*luhx\t  - h.)  ,6-2> 


and  so,  as  with  aiQ,  IiQ  is  a function  of  6^ 


ho  ■ <V, 


if 


K,  .1  o . n <1  . . V • r 

hsi 1 if  il  if 


*il]  \ 


■*)  ■ 


(6-13) 


With  the  definitions  of  equation  (6-8)  to  (6-11)  we  find  that  the  terms  involving 
the  deformation  modes  in  (6-7)  are 


XS"R  " E "i 


K - S7  A Q F 
♦if  ai0  ♦if 


- ST  PT  A 

♦if  6if  ai1 


(6-14) 
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The  other  terms  are  given  by  (6-5), 


■ E-. 


[ = > sT  I.nS.  - m.  < A2  + A A SA 

n L,  ^ *.f  10  *.f  x.f  x.f  *.f  a.Q 


+ S*  A A 

*if  ai0  *if  Xif 


(6-17) 


(6-18) 


7 THE  EQUATIONS  OF  MOTION 

Lagrange's  equations  merely  equate  the  generalised  forces  and  the 
generalised  effective  forces*.  Stating  this  equality  for  our  (n  + 6)  degrees  of 
freedom  under  the  assumption  of  small  perturbations  gives  two  sets  of  (n  + 6) 
simultaneous  equations.  With  the  datum  motion  that  we  have  taken  - briefly 
constant  linear  velocity,  zero  angular  velocity  in  a uniform  atmosphere 
(of  section  2)  - the  first  set  state  that  the  values  of  the  generalised  forces  in 
the  datum  motion  are  all  zero,  ie 


Qrf  = 0 


r = 1 -*■  (n  + 6)  (7-1) 


or  in  terms  of  the  separate  contributions  (aerodynamic,  gravitational,  etc) 


G^  + P.+E.  = 0 

rf  rf  rf 


r = 1 ->  n + 6 . (7-2) 


Expressions  for  the  individual  elements  in  these  equations  have  been  obtained 
above  - equations  (4-93),  (5-23),  (5-45)  and  (5-39). 

The  second  set  of  equations  express  the  continued  satisfaction  of 
D'Alembert's  principle  during  small  perturbations  of  the  datum  motion.  They  are 
written  as  the  matrix  equation 


* Often  they  are  expressed  on  the  equality  of  (on  the  left  hand  side)  the 
generalised  effective  forces  minus  the  conservative  generalised  forces,  and 
(on  the  right  hand  side)  the  non-conservative  generalised  forces.  This  is  the 
form  when  the  equations  are  written  in  terms  of  the  Lagrangian  function. 
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{[Ats]»2  * [C„]  * [Pj  * [E„]  ' [Qrs]} 


= o (7-3) 


where  the  constituent  matrices  are: 

(i)  the  inertia  matrix  [A  3»  equation  (6-7),  where  the  various  sub- 
matrices are  given  by  equations  (6-16),  (6-14),  (6-15),  (6-17)  and 
(6-18); 

(ii)  the  gravitational  matrix  [g  J,  equation  (5-25),  where  the  non-zero 
submatrices  are  given  by  equations  (5-26)  to  (5-28); 

(iii)  the  propulsive  matrix  [P^J,  equation  (5-45),  where  the  various  sub- 
matrices are  given  by  equations  (5-47)  to  (5-51); 

(iv)  the  structural  matrix  [e  ],  equation  (5-40),  which  has  one  non- 
zero submatrix  given  by  equation  (5-41); 

(v)  the  aerodynamic  matrix  - [Qrs]>  equations  (4-94)  and  (4-100), 
see  also  sections  4.5,  4.5.1  and  4.5.2. 

Equation  (7-3)  is  almost  the  same*  as  the  equation  given  in  Table  3 of 
Ref  1 except  that  it  has  been  written  in  terms  of  sectional  properties  and  in 
terms  of  modal  functions  appropriate  to  a sectional  description  of  the  configura- 
tion. In  particular  a minor  restriction  is  imposed  on  the  form  of  deformation 
in  that  what  might  be  called  ’chordwise  deformation'  is  forbidden  apart  from  that 
due  to  a flap  rotation;  but  more  noticeable  are  the  additional  terms  in  the  qq 
submatrices  resulting  from  different  ways  used  to  express  the  deformation 
(of  Appendix  A).  Thus  in  Ref  1 the  submatrix  G '■  , for  example,  is  null, 

whereas  in  the  present  development  it  is  given  by  equation  (5-26) . There  are 
consequently  some  corresponding  differences  in  the  transformation  to  other  forms 
of  the  equations  of  motion.  These  differences  have  been  demonstrated  in  detail 
in  Appendix  C. 


* When  the  choice  of  constant-velocity  axes  is  the  same.  In  Ref  1 they  have 
been  chosen  to  coincide  with  the  principal  axes  of  inertia  during  the  datum 

motion.  This  is  equivalent  to  putting  ^*m.A_  = 0 , and 


z_ , 


l X. 


ifO 


1^  = diagjl^  I I^{,  if  the  present  development. 
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8 CONCLUDING  REMARKS 

This  paper  has  been  written  as  a companion  document  to  Ref  1 . In  each  case 
the  equations  of  motion  of  an  aircraft,  for  small  perturbations  from  flight  with 
constant  linear  and  zero  angular  velocity,  have  been  developed  in  detail.  The 
constrast  has  been  that,  whereas  in  Ref  1 we  took  an  overall  view  of  the  aircraft, 
we  have  here  taken  a fore  and  aft  strip  of  the  aircraft  as  our  basic  unit,  con- 
sidered the  forces,  etc  on  the  strip,  and  built  up  from  that.  As  a consequence, 
the  representation  of  the  aircraft  deformation  in  the  present  paper  cannot  be 
made  to  correspond  exactly  with  that  of  the  earlier  paper.  The  deformat ional 
representation  is  basically  more  complicated;  and  other  complications,  such  as 
local  axes  for  each  strip,  are  also  introduced  by  strip  approach.  Of  course  a 
lot  of  these  complications  will  disappear  in  the  simplest  cases  (c/  Ref  5)  but 
even  so,  one  would  not  recommend  the  use  of  the  present  method  unless,  as  may 
quite  possibly  be  the  case,  it  means  that  adequate  basic  data  can  be  obtained 
much  more  simply,  as  for  example  if  two-dimensional  aerodynamic  theory  can  be 
used. 
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Appendix  A 

A GENERALISATION  OF  THE  GENERALISED  FORCE  EXPRESSION  OBTAINED  IN  REF  1 


In  Ref  1 it  was  assumed  (equation  (1))  that  the  deformation  of  the  air- 
craft was  given  precisely  by  a first  order  expression,  viz: 


[>>1 

n 

= 

xf 

+ R 

ql 

^n> 

yf 

• 

z<n) 

n 

_zf_ 

_qn 

This  assumption  simplified  some  of  the  subsequent  development,  while  making  no 
significant  restriction  to  the  freedom  to  represent,  to  a good  approximation, 
any  deformation  of  the  aircraft.  However,  in  certain  cases,  such  as  a deforma- 
tion which  involves  rotation  of  part  of  the  aircraft  as  a rigid  body,  it  is  not 
possible  to  represent  the  deformation  exactly  by  (A-l)  and  still  keep  the  number 
of  degrees  of  freedom,  n , finite.  The  deformations  assumed  in  this  paper  are 
indeed  of  this  type,  and  with  the  chosen  degrees  of  freedom,  their  representation 
is  (of  section  3) 
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(A- 2) 


If  we  make  use  of  the  expansions 
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along  with  the  relationships 
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on  the  flap  part  of  the  strip 
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This  equation  is  a refinement  of  ( 3— 18).  It  will  be  seen  that  the  additional 
(of  with  Ref  1)  second  order  term  has  the  form 
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Thus,  from  equation  (73)  of  Ref  1 the  matrix  of  generalised  coefficients 
resulting  from  a typical  local  force  vector 
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For  convenience  we  repeat,  from  Ref  1 (equation  (22))  the  formula  for  the  datum 
motion  values  of  the  generalised  forces,  viz: 
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(A- 13) 


The  additional  second  order  terms  (A-8)  in  the  expressions  for  the  deformation 
(A-7)  will  of  course,  make  no  difference  to  the  first  order  approximations  to 
the  overall  forces  and  moments  on  the  aircraft  and  so,  from  Ref  1 (equations 
(74)  to  (77)) 
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Appendix  B 

DETERMINATION  OF  THE  LOCATION  OF  THE  PRINCIPAL  AXES  OF  INERTIA 


The  determination  of  the  location  of  the  principal  axes  of  inertia  of  the 
aircraft  in  the  datum  state  is  basically  simple  but  for  completeness  it  was  con- 
sidered worthwhile  to  restate  it  here. 

Let  us  take  an  arbitrary  set  of  body-fixed  axes  as  a first  guess  at  the 
principal  axes  of  inertia.  Having  done  so  we  can  then  divide  the  aircraft  into 
strips,  locate  the  strip-fixed  axes  for  each  strip,  and  evaluate  the  quantities 
etc,  cj>^£  etc,  6^  for  each  strip  (af  section  3).  Note  that  under  our 
definition  the  strip-fixed  axes  are  found  by  moving,  in  one's  imagination,  the 
flap  part  of  the  strip  until  the  position  is  found  where  the  line  joining  the 
eg  of  the  strip  and  the  eg  of  the  flap  part  of  the  strip  passes  through  the 
hinge.  This  line  then  gives  the  direction  of  the  x strip-fixed  axis.  A point 
is  chosen  on  this  line  as  the  origin  of  these  axes  (ie  the  strip  reference  point), 
and  the  y axis  is  'drawn',  as  near  as  possible,  normal  to  both  the  9ides  of  the 
strip . 

One  can  then  evaluate  the  quantities  (af  equations  (6-5),  (6-6),  (6-18) 
and  (6-19)) 


xH  (say) 


■n 


6m  fx 


(B-l) 


Ib  (say)  = 


2 2 

6mA  + mA 

- Xf  *b 


(B-2) 


If  our  guess  of  axes  had  been  correct,  then  jx^  y^  z^f  would  be  zero  and  1^ 

would  be  diagonal.  If  this  is  not  so,  then  we  find  the  axes  transformation 
. T 

matrix  which  makes  diagonal.  Thus  the  rows  of  are  the 

eigenvectors  of  I normalised  so  that  . This  however  leaves  some 

ambiguity  - there  are  six  possible  matrices  - due  to  the  facts  that  any  of 

the  three  principal  axes  could  be  called  the  x-axis,  and  having  done  so,  there 
is  still  the  choice  to  be  made  as  to  which  is  the  positive  x direction.  One 
suggests  therefore  making  the  choice  which  makes  the  three  diagonal  elements  of 
S(jib  most  nearly  unity.  The  first  row  and  last  column  of  S$b  are  as  shown  below 
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cos  9.  cos  <p, 
b b 


cos  6.  sin  ip, 
b b 


- sin  0 

b 

cos  9,  sin  <p, 
b b 

cos  8,  cos  <i>, 
b b 


(B-3) 


and  so,  when  S^,  is  known,  one  can,  taking  the  principal  values,  in  general*, 

uniquely  determine  the  three  Euler  angles  <f>,  , 0,  , ip.  . The  translation  x.  , 

D D u b 

y^,  z^  , followed  by  the  Euler  rotations  0b>  » will  then  move  our  guessed 

axes  into  coincidence  with  the  principal  axes  of  inertia  of  the  aircraft.  Of 
course,  one  does  not  really  need  to  determine  the  angles  ^ etc;  the  matrix 
S<j,^  suffices  to  give  the  orientation  of  the  principal  axes  of  inertia. 

When  the  location  and  orientation  of  the  principal  axes  of  inertia  has 
been  determined  one  repeats  the  procedure: 

(i)  Divide  the  aircraft  into  strips. 

(ii)  Locate  the  strip-fixed  axes  for  each  strip. 

(iii)  Evaluate  for  each  strip 


if 

yif  Zif 

if 

6if  *if 

if 

* The  exception  is  when  cos  0^  = 0 . It  then  turns  out  that  one  can  get  only 
either  the  sum  or  the  difference  of  and  'I'b  • This  does  not  matter  of 

course  for  this  is  the  case  when  the  carried  axis  about  which  the  third  rota- 
tion is  made  coincides  with  the  original  axis  about  which  the  first  rotation 
was  made. 
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Appendix  C 

A BRIEF  CONSIDERATION  OF  THE  USE  OF  STRIP  THEORY 
IN  A BODY-FIXED  AXES  CONTEXT 


To  use  Lagrange's  equation  for  a non-inertial  frame  one  requires  that  the 
generalised  coordinates  for  which  it  is  used  must  not  influence  directly  the 
position,  etc  of  the  non-inertial  frame.  Taking  the  non-inertial  frame  to  be  the 
body-fixed  axes,  and  denoting  the  strip  in  which  the  body-fixed  axes  are  fixed 
as  strip  0,  we  therefore  choose  generalised  coordinates  -+■  for  the 
deformational  freedoms  such  that  the  aircraft's  perturbation  from  its  datum 
state  can  be  achieved  by  the  following  successive  steps: 

„ (c)  - (c)  • 

(i)  Translations  x^  , y^  , , in  the  directions  of  the  constant- 

velocity  axes  (also  called  9n+]  "*■  9n+3  respectively). 

(ii)  Euler  rotations  $,  §,  $ (in  that  order),  about  the  carried  (body- 
fixed)  axes  (also  called  9n+^>  qn+5’  qn+4  respectively). 

(iii)  Deformations  relative  to  the  body-fixed  axes  such  that  the  position 
and  orientation  of  the  strip-fixed  axes  of  strip  i , relative  to 
the  body-fixed  axes,  are  given  by  (of  equations  (3-2)  and  (3-3)) 


(C-1) 


In  addition,  the  flap  angle  perturbation  is  (of  equation  (3-7)) 


96 


Appendix  C 


Consequently,  due  to  a typical  force  distribution  which  is  such  that  the  overall 
forces,  moments  and  hinge  moment  on  the  ith  strip  are 
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*n+6 


the  contributions  to  the  generalised  forces  for  the  deformat ional  freedoms  are, 
calculated  by  the  principle  of  virtual  work  assuming  the  body-fixed  axes 
stationa ry , 
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Now  (of  equations  (3-13),  (3-15)  and  (4-85)) 
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(C-9) 


(C-10) 


(C-U) 


Lagrange's  equations  for  a non-inertial  frame,  taking  the  non-inertial  frame 
to  be  the  body-fixed  axes,  will  only  provide  us  with  the  n equations  for  the 
deformational  freedoms.  The  other  six  equations  are,  as  in  Ref  1,  obtained  from 
consideration  of  the  rates  of  change  of  the  linear  momentum,  and  angular 
momentum  about  the  aircraft  reference  point,  resolved  along  the  body-fixed  axes. 
Thus  we  require  the  resolutes  along  the  body-fixed  axes  of  the  total  applied 
force  and  moment  about  the  reference  point.  The  contributions  to  these  from  a 
typical  force  distribution  and  from  the  ith  strip  are  (of  equations  (C-4)  and 
(C-5)) 
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(C— 13) 


It  is  not  as  easy  to  relate  the  derivation  of  this  Appendix  to  that  for 
constant-velocity  axes  (in  the  main  part  of  this  paper)  as  it  was  in  Ref  1 when 
the  deformation  was  not  written  in  a form  attractive  for  strip  theory. 
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The  reader  will  have  noted  that  in  equations  (C— 1 ) and  (C-2),  we  could  not  put 
{k  - (K)q|  and  {f  - (F ) q J » where  K and  F represented  the  same  modal  functions 
as  in  the  main  part  of  the  paper,  but  had  to  add  the  distinguishing  circumflexes. 
In  particular  it  will  be  appreciated  that  the  term  (F)q  jcjj  ...  9n}»  in  (C-2), 
does  not  represent  a rigid  body  rotation,  not  even  to  a first  approximation. 

It  is  an  approximately  constant  rotation  about  axes  which  vary  from  strip  to 
strip,  or  in  other  words  it  is,  about  one  set  of  axes,  the  rotation 


if 


+ higher  order  terms 


(C— 14) 


which  varies  from  strip  to  strip  ( [of  equations  (3-13)  and  (3-15)).  From 
equation  (4-7)  it  will  be  seen  that  this  rotation  is  approximately  constant  if 
is  the  same  for  all  strips. 

If  the  deformation  of  section  3 is  associated  with  just  sufficient  rigid 
body  motion  to  make  the  displacement  and  slope,  at  the  aircraft  reference  point, 
due  to  the  deformation,  zero,  it  means  that  the  rigid  body  displacements  are: 


rotations  4>  satisfying* 
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and  translations* 
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* This  means  there  is  no  change  in  orientation  or  position  in  going  from  the 

unperturbed  strip-fixed  axes  to  the  strip-fixed  axes,  for  strip  0 (af  equations 
(4-38)  and  (4-39)). 
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where  P 

q 

Thus 


is  given  by  equation  (5-55). 
if  we  take 
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K - 


which  for  arbitrary  F , K can  be  achieved  by  taking 


F = F - Q"1  S.  Si  Q.  <F)n 
*if  *if  *0f  *0f  ° 


(and  so,  as  expected,  taking  (F)Q  = 0),  and 


K = K + A P 
Xif  * 


(C-20) 


(C— 2 1 ) 


(C-22) 


then,  to  a first  approximation,  the  modes  of  deformation  relative  to  the  strip- 
fixed  axes  will  be  the  same  with  both  representations  (that  of  section  3 and  that 
cf  this  Appendix).  This  suggests  that  with  the  modal  relationships  given  above 
the  following  relationship  might  not  be  far  from  the  truth 

■ T1  -(«)«+  A.  pJT  . (023) 
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From  (5-11)  and  (C— 11)  it  is  easily  seen  to  be  true  for  the  datum  motion  values  of 
these  generalised  forces.  For  both  models  to  represent,  to  a good  approximation, 
the  same  perturbation  there  will  be  a similar  relationship,  corresponding  to 
(C-23),  between  the  generalised  coordinates,  which  we  can  deduce  to  be 


0 0 q 


1®0*V«I  1 0 

-p  oi 


(C-24) 


Comparing  the  expressions  for  the  overall  forces  moments  and  hinge  moment  on  a 
strip  - equations  (C-4)  to  (C-6)  with  equations  (5-7)  to  (5-9)  we  then 
immediately  obtain  the  relationships 


X..  ...  X.  , 
ll  i,n+6 


X.,  ...  X. 
ll  i,n+6 


- {(K)0 


+ A P 

x0f  qI 


(C-25) 


when  the  modal  functions  are  related  by  equations  (C— 19)  and  (C-20) . 

ing  these  relationships  into  equations  (C— 12)  and  (C-13),  and  writing 


* <w  + 


11L?  r~ 


n+h , fj 


k qk 


(C-26) 


n+6,k 


we  then  find,  after  some  analysis,  that  (.of  equations  (5-11)) 
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^n+1 ,f 

s 

Q ( i ) 

Vt-l.f 

• 

_ 

• 

• 

Q(i) 

Jn+6,f_ 

^n+6,f_ 

(C-27 ) 


and 


"o(i) 

6(i) 

Vi.i 

,n+6 

^n*6, 1 

_ 

o i o 

0 0 I 


- ®.  * V'l 


- P 


0 0 

I 0 

0 I 


SI  A /-  (oOft*A  P I 

•if  x!“s)  ^if  l ° x0f  ql 

if 


♦ X 


A ST 

’x<“,T 

♦ ST 

V::,r 

‘if  *if 

if 

:(us) 

if 

j(if) 

♦if 

if 

M(U,) 

Mif 

ij<Ui) 

"if 

L if  - 

A S* 
\f  *il 


J<u*^ 

if 

;(u«) 

if 

;<ut)| 

.if 


r-(u.)|\ 

Lif 

=(u»)i 

Mif 
5(»»>i 
.if 


(C-28) 


This  is  exactly  the  same  relationship  as  that  obtained  in  Ref  1 (equation  (161)). 
The  relationship  of  Ref  1 also  holds  for  those  coefficients  which  express  the 
influence  of  the  body  freedoms  on  the  deformational  freedoms,  but  not,  as  one 
would  expect  because  of  the  differences  displayed  in  Appendix  A,  for  the 
'deformation-deformation'  coefficients.  Comparing  (C-ll)  and  (5-12)  we  find  that, 
for  the  modal  relationships  of  (C— 19)  and  (C-20), 
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In  Ref  1,  because  of  the  slightly  different  model,  the  expression  for  & is 
different  and  is  & = where  (of  equation  (161)  of  Ref  1) 


& 


77 


PTST 

qSi 


if  x 


(us)' 


if 


II' 


+ > ( rV  A . .S,  P - PTS^  A . ,S,  R 

*if  e<“s)  ♦if  * ♦if  i<“s)  ♦if 

if  if 


PTS^  A . .S 
q *.f  ~(us)  * 

if 


(k  + sj  A , .S.  V 1 

l,Xif  ♦if  xs(->  ♦ifj 


(c- 31) 


the  summation  is  over  the  ith  strip,  and  R is  given  by  equation  (3-20).  One 
term  is  clearly  common  to  (C-30)  and  (C-31). 
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GLOSSARY  OF  TERMS 


(i)  Frames  of  reference  (all  right-handed  orthogonal  cartesian) 


body-fixed  axes 


axes  whose  origin  and  orientation  are  fixed  in  a 
small  material  portion  of  the  body.  (The  ones  used 
are  such  that  the  origin  is  in  the  main  part  of 
strip  0 during  the  datum  motion.)* 


constant-velocity 


axes  having  constant  linear  and  angular  velocity 
relative  to  an  inertial  frame.  (The  ones  used  have 
zero  angular  velocity  and  are  taken  to  be  coincident 
with  the  datum-attitude  earth  axes  during  the  datum 
motion. ) 


datum-attitude 


axes  with  which  the  body-fixed  axes  coincide  during 


earth  axes 


the  datum  motion 


datum-path  earth 


axes  with  the  x-axis  in  the  direction  of  the  datum 


velocity,  the  xz  plane  vertical,  and  the  z axis 


downwards 


no-deformation-body- 


fixed  axes 


an  axes  system,  arbitrary  except  insofar  as  it  is  of 
the  same  order  of  nearness  to  the  body-fixed  axes 
and  the  datum-attitude  earth  axes.  (These  three 


frames  of  reference  are  all  assumed  to  remain  close 


together  during  the  perturbed  motion.) 


normal  earth-fixed 


axes  fixed  relative  to  the  earth  with  the  z-axis 


vertically  downwards 


principal  axes  of 


inert  la 


axes  with  origin  at  the  centre  of  gravity  of  the  body 
and  such  that  the  three  products  of  inertia  about 


the  axes  are  zero 


strip-fixed  axes 


axes  fixed  in  a strip  of  the  aircraft  such  that  their 
origin  is  at  the  strip  reference  point  the  x-axis 
passes  through  a point  on  the  strip  hinge  (and  in  the 
unperturbed  state,  the  strip  eg),  the  y-axis  is 
approximately  normal  to  the  planar  sides  of  the  strip, 


* Whether  the  body-fixed  axes  remain  orthogonal,  during  perturbations  of  the 
aatum  motion,  or  not,  is  irrelevant  to  the  present  development. 
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GLOSSARY  OF  TERMS  (continued) 


strip-fixed  axes  and  the  positive  x-direction  is  towards  the  strip 

(continued)  leading  edge  - of  Appendix  B.  (For  a strip  without 

a hinge  one  would  take  the  x-axis  to  be  some  chord 

line  and  the  y-axis  to  be,  as  near  as  possible, 
normal  to  the  planes  dividing  the  strip  from  the 
adjacent  strips.) 

unperturbed-strip-  axes  whose  position  coincides  with  the  position 

fixed  axes  the  strip-fixed  axes  would  occupy  if  there  was  no 

perturbation  from  the  datum  motion. 

(ii)  Orientation 

(a)  Attitude  angles  defining  the  attitude  of  the  aircraft  relative  to  the 
normal  earth-fixed  axes. 

angle  of  bank  angle  between  the  z-axis,  of  the  body-fixed  axes  and 

the  vertical  plane  containing  the  x-axis  of  the 
same  frame 

angle  of  inclination  angle  between  the  x-axis  of  the  body-fixed  axes 

and  the  horizontal  plane 

nose-azimuth  angle  angle  between  the  projection  of  the  x-axis,  of  the 

body-fixed  axes,  on  the  horizontal  plane,  and  the 
x-axis  of  the  normal  earth-fixed  axes. 

(b)  Flight  path  angle  defining  the  direction  of  flight  relative  to  the 
normal  earth-fixed  axes. 

angle  of  climb  angle  between  the  tangent  to  the  flight  path  and 

the  horizontal  plane,  positive  when  climbing 

angle  of  track  angle  between  the  x-axis  of  the  normal  earth-fixed 

axes  and  the  projection  on  the  horizontal  plane  of  the 
tangent  to  the  flight  path,  positive  when  a right- 
handed  rotation  about  the  downward  vertical  takes 
one  from  the  normal  earth-fixed  x-axis  to  the  flight 
path  projection. 
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GLOSSARY  OF  TERMS  (concluded) 

(c)  Incidence  angles  defining  the  direction  of  flight  relative  to  the  air 


incidence  magnitude*  magnitude  of  angle  between  the  x-axis  of  the  body- 

fixed  axes  and  the  direction  of  the  velocity  of  the 
aircraft  relative  to  the  air 


incidence-plane  angle  between  the  incidence  plane  and  the  zx  plane 

angle*  of  the  body-fixed  axes. 


(iii)  Miscellaneous 

aircraft  reference  the  material  point  on  the  aircraft  which  is  the 
point  origin  of  the  body-fixed  axes 


attitude  deviation  the  matrix  which  by  premultiplying  a column  vector 
matrix  changes  the  axes  directions  in  respect  to  which  the 

vector  is  resolved  into  components 


axes  transformation  synonymous  with  attitude  deviation  matrix 
matrix 


incidence  plane 


reference  axis 


reference  point 

strip  reference 
point 


the  plane  defined  by  the  x-axis  of  the  body-fixed 
axes  and  the  direction  of  the  velocity  of  the  aircraft 
relative  to  the  air 

a line,  not  necessarily  straight,  joining  the  strip 
reference  points  of  the  strips  of  a component,  such 
as  the  wing,  of  the  aircraft 

synonymous  with  aircraft  reference  point 

a chosen  material  point  on  the  main  part  of  a strip 
taken  as  the  origin  of  the  strip-fixed  axes. 


* These  definitions  can  similarly  be  used  for  a strip  or  aerofoil  section, 
in  terms  of  the  strip-fixed  axes. 
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L'lST  of  symbols 


A,  etc 


A., 

jk 


B • • , B . ; etc 

IX  10 

B.  . ,B. • . etc 

1XJ  1XJ 

B* 

i0J 

B. . ,B. .Q  etc 
ljx’  1J0 


- E. 
J 

E., 

Jk 


- G. 
J 


Jk 

H 

I 

I ,1  ,1 

x y z 


Ii0 


Xil 


J,  etc 
<P 

K 

L 

L 

L , L 
8 P 

L.,L  . etc 
i gi 

M 

M 


skew-symmetric  matrices  involving  <{>,  0,  <P,  etc  (see 
equation  ( 3— 16)) 

ponderous  inertia  coefficient 

aerodynamic  hinge  moment  on  the  ith  strip 

aerodynamic  hinge  moment  coefficients  for  ith  strip 

aerodynamic  hinge  moment  coefficients  of,  say,  a velocity  of 
strip  j , in  the  expression  for  the  hinge  moment  on  strip  i 
(2) 

a modified  B£qj  “ see  equation  (4-146) 

aerodynamic  hinge  moment  coefficients  expressing  the  effect,  on 
the  hinge  moment  on  strip  i , of  a change  in  relative  position 
and  orientation  of  strips  i and  j 

differential  operator  d/dt 

generalised  structural  force 

structural  stiffness  coefficient 

torsional  modal  matrix  (see  equation  (3-3)) 
generalised  gravitational  force 

gravitational  stiffness  coefficient 

Heaviside  step  function 
unit  matrix 

principal  moments  of  inertia  of  undeformed  aircraft 

matrix  of  moments  and  products  of  inertia  of  undeformed  aircraft 
(see  equation  (6-6)) 

matrix  of  moments  and  products  of  inertia  of  the  ith  strip 
(see  equation  (6-8)) 

matrix  of  moments  and  products  of  inertia  of  the  flap  of  the 
ith  strip  (see  equation  (6-9)) 

matrices  formed  from  the  elements  of  |<j>  0 f etc, 
see  equation  (4-88) 

flexural  modal  matrix  (see  equation  (3-2)) 
aerodynamic  rolling  moment 
typical  rolling  moment 

gravitational,  propulsive  rolling  moments 

rolling  moments  on  ith  strip 

aerodynamic  pitching  moment 
typical  pitching  moment 
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Mq.M, 

M ,M 
g P 

M.,M  . etc 
i gi 

M. • ,M. x etc 

IX  10 

M.  • . ,M. • . etc 
lxj  iz j 

M*  . ,M*  . 
lxj’  1ZJ 

M. . ,M. .Q  etc 

1JX*  1J0 


N 

N 

VNi 

N ,N 
g P 

N.,N  . etc 
i gi 

- P. 


jk 


Qi 


,(i) 


J 

Qjk 


R. 
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LIST  OF  SYMBOLS  (continued) 
modal  matrices,  see  Appendix  A 
gravitational,  propulsive  pitching  moments 
pitching  moments  on  ith  strip 

aerodynamic  pitching  moment  coefficients  for  ith  strip 

aerodynamic  pitching  moment  coefficients  of,  say,  a velocity  of 
strip  j , in  the  expression  for  the  pitching  moment  on  strip  i 

modified  M^xj  etc  ~ see  equations  (4-140)  to  (4-142) 

aerodynamic  pitching  moment  coefficients  expressing  the  effect, 
on  the  pitching  moment  on  strip  i , of  a change  in  the  relative 
position  and  orientation  of  strips  i and  j 

aerodynamic  yawing  moment 

typical  yawing  moment 

modal  matrices,  see  Appendix  A 

gravitational,  propulsive  yawing  moments 

yawing  moments  on  ith  strip 

generalised  propulsive  force 

matrix  of  modal  slopes  at  reference  point  (see  equations  (5-54) 
and  (5-55)) 

matrix  which  is  '0  factor'  of  axes  transformation  matrix  (see 
equation  (3-11)) 

propulsive  stiffness  coefficient 
generalised  aerodynamic  force 

total  generalised  force,  or  typical  contribution  to  generalised 
force 

matrix  relating  angular  velocities  and  orientation  (see 
equations  (3-14)  and  (4-3)) 

contribution  to  generalised  aerodynamic  force  from  strip  i 

typical  contribution  to  generalised  force  from  strip  i 

aerodynamic  coefficient 

modal  matrix  (see  equation  (3-20)) 

matrix  which  is  '<)>  factor'  of  axes  transformation  matrix  (see 
equation  (3-10)) 


1 10 


LIST  OF  SYMBOLS  (continued) 


U 

V 

W 

X 

X 

X ,x 
g p 

X.,X  . etc 
i gi 

X.*,X.jt  etc 

IX  10 

X.  • . ,X. • . etc 
ixj’  1ZJ 

X* 

X.  . ,X. ..  etc 
ijx’  1J0 


Y 

Y 

Y ,Y 

g P 


Y .  , Y . etc 
i gi 

Z 

Z 

Z ,Z 
g P 

Z. ,Z  . etc 
i gi 

Z.»,Z.;  etc 

IX  10 

Z . • . ,Z . • . etc 

ixj  1ZJ 


Zi0j 


axes  transformation  matrix  (see  section  3).  The  subscript  is 
added  when  it  is  necessary  to  specify  the  arguments  - in  this 
case  <Jk,  0^,  iJk  . Absence  of  a subscript  means  the  arguments 

are  <p,  0,  ip 

modal  matrix  (see  Appendix  A) 
airspeed 

kinetic  energy  relative  to  the  constant-velocity  axes 
overall  aerodynamic  force  resolute 
typical  overall  force  resolute 

gravitational,  propulsive  overall  force  resolutes 
overall  force  resoluted  on  strip  i 

aerodynamic  force  resolute  coefficients  for  ith  strip 

aerodynamic  force  resolute  coefficients  of,  say,  a velocity  of 
strip  j , in  the  expression  for  the  overall  force  on  strip  i 

modified  “ see  equations  (4-128)  and  (4-129) 

aerodynamic  force  resolute  coefficients  expressing  the  effect, 
on  the  overall  force  on  strip  i , of  a change  in  the  relative 
position  and  orientation  of  strips  i and  j 

overall  aerodynamic  force  resolute 

typical  overall  force  resolute 

gravitational,  propulsive  overall  force  resolutes 

matrix  which  is  '<(>  factor'  of  axes  transformation  matrix  (see 
equation  (3-12)) 

overall  force  resolutes  on  strip  i 

overall  aerodynamic  force  resolute 
typical  overall  force  resolute 

gravitational,  propulsive  overall  force  resolutes 
overall  force  resolutes  on  strip  i 

aerodynamic  force  resolute  coefficients  for  the  ith  strip 

aerodynamic  force  resolute  coefficients  of,  say,  a velocity  of 
strip  j , in  the  expression  for  the  overall  force  on  strip  i 

modified  “ see  equations  (4-130)  and  (4-131) 
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LIST  OF  SYMBOLS  (continued) 


I 


3 , . y d • i y3i « 

lCr  il*  i3 


ei* 


e . 
gi 


f . 
gi 


8i’V 


l , l ,Z‘ 

a 9 z * z 


m 


m. 

l 

6m 

n 

P 


q 

qi 

qr 


r 


r . 

l 


t 

u 


u . 

l 

u . 
mi 

v 


v. 

l 

v . 
mi 

w 


w. 

l 


w . 


mi 


certain  column  vectors  - see  equations  (5-19)  to  (5-21) 

x-component  of  local  aerodynamic,  gravitational  force  vector 
at  ith  strip 

flap  modal  vector  (see  equation  (3-7)) 

y-component  of  local  aerodynamic,  gravitational  force  vector 
at  ith  strip 

acceleration  due  to  gravity 

z-component  of  local  aerodynamic,  gravitational  force  vector 
at  ith  strip 

third  column  of  (see  equation  (5-14)) 

aerodynamic  derivatives  - obsolescent  notation  (see  section 
4.22) 

mass  of  aircraft 

mass  of  strip  i of  the  aircraft 
mass  of  a particle 

number  of  deformat ional  degrees  of  freedom 

angular  velocity  resolute 

angular  velocity  resolute  of  ith  strip 

angular  velocity  resolute 

angular  velocity  resolute  of  ith  strip 

generalised  coordinate 

angular  velocity  resolute 

angular  velocity  resolute  of  ith  strip 

time 

linear  velocity  resolute 

linear  velocity  of  ith  strip  reference  point 
particle  velocity  resolute  at  ith  strip 
linear  velocity  resolute 

linear  velocity  resolute  of  ith  strip  reference  point 
particle  velocity  resolute  at  ith  strip 
linear  velocity  resolute 

linear  velocity  resolute  of  ith  strip  reference  point 
particle  velocity  resolute  at  ith  strip 


I 


‘ 1 


J 
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LIST  OF  SYMBOLS  (continued) 


X1  (=  Vl} 


~i2 

^si 


particle  position  resolute 

resolute  of  reference  point  position,  relative  to  its  position 
in  datum  motion,  excluding  deformational  contribution 

resolute  of  strip  reference  point  position 

x-coordinate  of  strip  eg,  referred  to  strip-fixed  axes,  when 
the  flap  angle  is  zero 

x-coordinate  of  strip  eg,  in  datum  motion,  in  body-fixed  axes 
frame  of  reference 

x-coordinate  of  strip  hinge,  referred  to  strip-fixed  axes 


see  equation  (4-109) 


yl  (=  V2} 


particle  position  resolute 

resolute  of  reference  point  position,  relative  to  its  position 
in  datum  motion,  excluding  the  deformational  contribution 

resolute  of  strip  reference  point  position 


ifO 


y-coordinate  of  strip  eg,  in  datum  motion,  in  body-fixed  axes 
frame  of  reference 


Z1  (=  V3} 


'ifO 


particle  position  resolute 

resolute  of  reference  point  position,  relative  to  its  position 
in  datum  motion,  excluding  the  deformational  contribution 

z-coordinate  of  strip  eg,  in  datum  motion,  in  body-fixed  axes 
frame  of  reference 


see  equation  (4-110) 


E submatrix  of  structural  stiffness  coefficients 

W 

G ,G  .,G.  ,GX.  submatrices  of  gravitational  stiffness  coefficients 

qq  <?<)>  $q'  M 


L fL  ,L  , 
aq  ox  c< )> 


submatrices  of  coefficients  in  expression  for  aerodynamic 
rolling  moment  about  origin  of  constant-velocity  axes 


L ,L  ,L  . 
opq  cpx ’ op<p 


submatrices  of  coefficients  in  expression  for  propulsive  rolling 
moment  about  origin  of  constant-velocity  axes 


M.  ,M.  ,M.. 
tq  vx’  7<P 


submatrices  of  coefficients  in  expression  for  aerodynamic 
pitching  moment  on  ith  strip 


aq  * ox  e<Ji 


submatrices  of  coefficients  in  expression  for  aerodynamic  pitch- 
ing moment  about  origin  of  constant-velocity  axes 


M ,M  ,M  . 
opq  apx  opty 


submatrices  of  coefficients  in  expression  for  propulsive 
pitching  moment  about  origin  of  constant-velocity  axes 


N ,N  ,N  . 
oq  ox  otp 


submatrices  of  coefficients  in  expression  for  aerodynamic  yawing 
moment  about  origin  of  constant-velocity  axes 
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N ,N  ,N  x 
apq  opr  cp<p 

P »p  . P : 
qq  xq  pq 

P P 
x<p  ’ <t><t> 


submatrices  of  coefficients  in  expression  for  propulsive  yawing 
moment  about  origin  of  constant-velocity  axes 

submatrices  of  propulsive  stiffness  coefficients 


X ,X  ,X 
q ’ x*  P 

X.  ,X.  ,X.. 
tq  ix  ’ ^<f> 

X ,X  ,X  . 
pq  px’  pp 

q x’  $ 

y y y 

pq*  px*  pp 

Z ,Z  ,Z. 
q’  x*  4> 

Z.  ,Z.  ,Z.. 

tx’  t<)> 


Z ,Z  ,Z  . 
pq  px  pP 


submatrices  of  coefficients  in  expression  for  overall  aero- 
dynamic force  resolute 

submatrices  of  coefficients  in  expression  for  overall  aero- 
dynamic force  resolute  on  ith  strip 

submatrices  of  coefficients  in  expression  for  overall  propulsive 
force  resolute 

submatrices  of  coefficients  in  expression  for  overall  aerodynamic 
force  resolute 

submatrices  of  coefficients  in  expression  for  overall  propulsive 
force  resolute 

submatrices  of  coefficients  in  expression  for  overall  aero- 
dynamic force  resolute 

submatrices  of  coefficients  in  expression  for  overall  aero- 
dynamic force  resolute  on  ith  strip 

submatrices  of  coefficients  in  expression  for  overall  propulsive 
force  resolute 


r ,r  ,r  ,r, 

u w q f 


Ax,Az,A9 

© 

<1> 

* 


vortex  strength 

coefficients  in  expression  for  vortex  strength  (of  equation 
(4-13)) 

see  equations  (4-106)  to  (4-108) 
angle  of  inclination 
angle  of  bank 
nose-azimuth  angle 


a 


a 


0 


a 


l 


a . 
l 


0i 


incidence  of  aerofoil 
no  lift  incidence  of  aerofoil 

no  pitching  moment  incidence  for  aerofoil  - pitching  moment  is 
about  centre  of  the  circle  from  which  the  profile  is  generated 
by  conformal  transformation 

first  element  of  column  vector  given  by  equation  (3-13) 
second  element  of  column  vector  given  by  equation  (3-13) 
angle  of  climb 


Y 
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LIST  OF  SYMBOLS  (continued) 


third  element  of  column  vector  given  by  equation  (3— 13) 


6. 

l 


angle  of  flap  rotation  at  strip  i 


6m 


mass  of  a particle 


9 V5> 

e. 


e . 

U1 


orientation  angle  of  no-deformation-body-fixed  axes  relative 
to  constant-velocity  axes 

orientation  angle  of  strip-fixed  axes  relative  to  no-deformation- 
body-fixed  axes 

orientation  angle  of  strip-fixed  axes  relative  to  unperturbed- 
strip-fixed  axes 


see  equation  (4-111) 


P 

T 


^ (-  V4> 


air  density 

see  equation  (4-87) 

orientation  angle  of  no-deformation-body-fixed  axes  relative 
to  constant-velocity  axes 


orientation  angle  of  strip-fixed  axes  relative  to  no-deformation- 
body-fixed  axes 

orientation  angle  of  strip-fixed  axes  relative  to  unperturbed- 
strip-fixed  axes 


xif ,Xiu,Xiw 
Xiq’Xi6 

* W 


coefficients  in  expression  for  local  aerodynamic  force 
(of  equation  (4-15)) 


orientation  angle  of  no-deformation-body-fixed  axes  relative  to 
constant-velocity  axes 

orientation  angle  of  strip-fixed  axes  relative  to  no-de format ion- 
body-  fixed  axes 

orientation  angle  of  strip-fixed  axes  relative  to  unperturbed- 
strip-fixed  axes 


si 


e 

etc 


see  equation  (C-30) 
see  equation  (C— 3 1 ) 


operator  defined  by  equation  (4-127) 

operator  signifying  the  'steady  part  of’  (of  equation  (4-126)) 

lower  triangular  matrix  formed  from  the  elements  of  6 ty}  etc  - 
of  equation  (4-95) 
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LIST  OF  SYMBOLS  (continued) 


Dressings 


(i)  Subscripts 


A quantity  is  relative  to  or  about  the  origin  of  a particular  set  of 


axes  where 


c denotes  the  constant-velocity  axes 

n denotes  the  no-deformation-body-f ixed  axes 

s denotes  the  strip-fixed  axes  (this  subscript  is  often  omitted  to  avoid 
confusion  with  s for  structural) 

u denotes  the  unperturbed-strip-fixed  axes 

This  subscript,  when  present,  is  always  placed  first. 

i indicates  point  or  force  on  ith  strip 

g indicates  gravitational 

p indicates  propulsive 

s indicates  structural 

0 indicates  the  strip  i = 0 (the  one  containing  the  aircraft  reference 

point),  or  a leading  edge  suctional  force,  or  a quantity  associated 
with  the  strip  eg 

The  following  two,  when  present,  are  always  placed  last  except  when 
followed  by  a nought. 

e datum  value 

f value  during  datum  motion 


(ii)  Superscripts 
' applied  to 


applied  to 


indicates  a summation  over  the  main  part  of  a strip 
indicates  a summation  over  the  flap  part  of  a strip 


* denotes  certain  modified  coefficients 

T denotes  the  transpose  of  a matrix 

(2)  denotes  the  two-dimensional  value 
(i)  denotes  the  contribution  from  strip  i 

Other  bracketed  superscripts  denote  the  axes  of  resolution,  viz: 

(c)  constant-velocity  axes 

(dp)  datum-path  earth  axes 

(is)  strip-fixed  axes  of  the  ith  strip 


(dot)  denotes  derivative  with  respect  to  time 

(bar)  denotes  typical  or  total 

(circumflex)  refers  to  body-fixed  axes,  encastre  modes,  displacement 
body  freedoms 

(tilde)  see  equation  (A-8) 


TM  Str  933 


TM  Str  933 


No. 


Author 


1 D.L.  Woodcock 


2 D.L.  Woodcock 


3 H.R.  Hopkin 


4 D.L.  Woodcock 


5 J.C.A.  Baldock 


6 D.L.  Woodcock 


7 D.L.  Woodcock 


8 D.E.  Davies 


9 D.E.  Davies 


117 


REFERENCES 


Title,  etc 


Divers  forms  and  derivations  of  the  equations  of  motion 
of  deformable  aircraft  and  their  mutual  relations  hip. 
RAE  Technical  Report  77077  (1977) 

Several  formulations  of  the  equations  of  motion  of  an 
elastic  aircraft  as  illustrated  by  a simple  example. 

RAE  Technical  Memorandum  Structures  912  (1977) 

A scheme  of  notation  and  nomenclature  for  aircraft 
dynamics  and  associated  aerodynamics. 

ARC  R&M  3562  (1966) 

Mathematical  approaches  to  the  dynamics  of  deformable 
aircraft  - Part  II  The  dynamics  of  deformable  aircraft. 
ARC  R&M  3776  (1971) 

The  equations  of  motion  of  a flexible  aircraft 
symmetric  motion,  small  perturbations  from  trimned 
straight  and  level  flight. 

Unpublished  note  (1977) 

Formulation  of  the  equations  of  motion  of  a deformable 
aircraft  using  Lagrange's  equations  in  an  arbitrary 
non-inertial  frame  of  reference. 

RAE  paper  to  be  issued 

A formulation  of  the  equations  of  motion  of  a semi-rigid 
deformable  aircraft  when  only  the  deformations  are  small. 
RAE  Technical  Memorandum  Structures  914  (1977) 

Calculation  of  generalised  airforces  on  two  parallel 
lifting  surfaces  oscillating  harmonically  in  subsonic 
flow. 

ARC  R&M  3749  (1974) 

Generalised  aerodynamic  forces  on  a T-tail  oscillating 
harmonically  in  subsonic  flow. 

ARC  R&M  3422  (1966) 


Fig  2 Two-dimensional  tandem 


|/0/bS J 


1 


REPORT  DOCUMENTATION  PAGE 

Overall  security  classification  of  this  page 


UNLIMITED 


As  far  as  possible  this  page  should  contain  only  unclassified  information.  If  it  is  necessary  to  enter  classified  information,  the  box 
above  must  be  marked  to  indicate  the  classification,  e.g.  Restricted,  Confidential  or  Secret. 


1 . DR1C  Reference  2.  Originator’s  Reference  3.  Agency  4.  Report  Security  Classification/Marking 

(to  be  added  by  DRIC)  Tv,  Reference 

RAE  „ Q,,  H/.  UNLIMITED 

Structures  933  N/A 


5.  DRIC  Code  for  Originator 
850100 


6.  Originator  (Corporate  Author)  Name  and  Location 

Royal  Aircraft  Establishment,  Farnborough,  Hants,  UK 


5a.  Sponsoring  Agency’s  Code  6a.  Sponsoring  Agency  (Contract  Authority)  Name  and  Location 


7.  Title  The  use  of  strip  theory  in  the  dynamics  of  deformable  aircraft 


7a.  (For  Translations)  Title  in  Foreign  Language 


7b.  (For  Conference  Papers)  Title,  Place  and  Date  of  Conference 


8.  Author  1.  Surname,  Initials  9a.  Author  2 


Woodcock,  D.L. 


1 1 . Contract  Number 


12.  Period 


9b.  Authors  3,  4 .... 


13.  Project 


10.  Date  Pages  Refs. 

i "s  I » 


14.  Other  Reference  Nos. 


15.  Distribution  statement 

(a)  Controlled  by  - 

(b)  Special  limitations  (if  any)  - 


16.  Descriptors  (Keywords)  (Descriptors  marked  * are  selected  from  TEST) 

Aeroelasticity*.  Equations  of  motion*.  Dynamics*. 


17.  Abstract 

A detailed  formulation  of  the  equations  of  motion  of  a deformable  aircraft  is 
given.  The  development  is  from  Lagrange's  equations  for  an  inertial  frame,  and  is 
made  in  terms  of  the  position,  orientation,  force  and  inertia  properties  of  narrow 
strips  of  the  aircraft  which  lie  fore  and  aft  in  the  unperturbed  state.  The  latter 
is  one  of  constant  linear  velocity  and  zero  angular  velocity.  Particular  account  is 
taken  of  the  deformation  and  loading  in  the  unperturbed  state. 


RAE  Form  A 143 


